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ABSTRACT 



Aims. We aim both to complement the existing data on the spin history of the peculiar accreting X-ray pulsar GX 1+4 with more past 
and current data from BeppoSAX, INTEGRAL, and Fermi and to interpret the evolution in the framework of accretion theory. 
Methods. We used source light curves obtained from BeppoSAX/WFC and INTEGRAL/IS GRI to derive pulse periods using an 
epoch-folding analysis. Fermi/GBM data were analysed by fitting a constant plus a Fourier expansion to background- subtracted rates, 
and maximizing the Y 2 statistic. We completed the sample with hard X-ray light curves from Swift/BAT. The data were checked for 
correlations between flux and changes of the pulsar spin on different timescales. 

Results. The spin-down of the pulsar continues with a constant change in frequency, i.e., an apparently accelerating change in the 
period. Over the past three decades, the pulse period has increased by about ~50%. Short-term fluctuations on top of this long-term 
trend do show anti-correlation with the source flux. Possible explanations of the observed long-term frequency and its dependence on 
flux are discussed. 

Key words, accretion, accretion discs - X-rays: binaries - stars: pulsars: individual (GX 1+4) 



1. Introduction 

Accreting X-ray pulsars are highly magnetized neutron stars in a 
binary system, which accrete matter from their companion star. 
The mass transfer can take place via Roche-Lobe overflow for 
low mass X-ray binaries (LMXBs), strong stellar winds for giant 
stars in high mass X-ray binaries (HMXBs), or the Be emission 
mechanism for Be X-ray binaries. These accreting pulsars radi- 
ate predominantly in the X-ray band, and the radiation is mod- 
ulated by the st ellar rotation of the pulsar. For a review of this 
subject, see e.g.. lNagasel(ll989l) . 

Since the discovery o f the first a ccreting X-ray puls ars, Cen X-3 
dGiacconi et al.lll97ll). GX 1+4 (ILewin et al.ll 19711) . and Her X- 
1 ( Tananbaum et al.ll 9721) . the pulse periods of well-established 
accreting X-ray pulsars have been monitored more or less reg- 
ularly with a wide variety of high-energy observatories. Early 
observations mainly found sources being spun up, which was 
easily explained by the acceleration through the angular momen- 
tum of accreted matter. Further investigations have demonstrated 
a wide variety in the pulse-period evolution, with some sources 
best described by a random walk , others showing clear secu- 
lar changes. Various model s (e.g.. lGhosh & Lamblll977t IWand 
Il987h iLovelace et al.l[l9 95) have been proposed to explain the 
observed behaviour. 

GX 1+4 was discovered in 1970 by a balloon X-ray ob- 
servation at energies above 15 k eV showing pulsat ions with 
a period of about two minutes (ILewin et al.l 1 1 97 lb . During 
the 1970s, it was one of the brightest X-ray sources in the 



Gala ctic centre region and was globally spinning up strongly 
(e.g.. iDotv et al.l[l98Tl IWarwick et al.l [1981 IWhite et all [19831: 
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lMcChntock&Leventhallll989h . In the early 1980s the source 
went through a low state in X-ray flux and remained unde- 
tectable, because at least two orders of magnitude below the pre - 
viously observed levels (lHall & Davelaarlll983l: lMukailll988l) . 
When th e source was detected again it had undergone a torque 
reversal dMakishima et alJI 19881) . and ever since it has been gen- 
erally spinning down strongly, see AppendixlAl 

The optical counte rpart of GX 1+4 was discovered by 
lGlass&Feas1(fl97l as the infrared source V21 16 Oph. The op- 
tical composite emission spectrum of the proposed counterpart 
indicated that the object was almost certainly a binary system, 
consisting of a symbi ot ic red giant and a mu c h hott er source 
dDavidsen et al.l Il977l) . IChakrabartv & Rochd (1 19971) confirm 
the optical companion to be V21 16 Oph. More recently, infrared 
observations have indicated a mass of about 1 .2 M for the M gi- 
ant star (assuming a mass of about 1.35 M for the neutron star), 
implying that the M giant star is a first ascent giant that does not 
fill its Roche Lobe dHinkle et alJl2006l) . Therefore, GX 1+4 is a 
LM XB that is capturing the stellar wind of its M6 III compan- 
ion ( IChakrabartv & Rochd[T997l:lHinkle et al.ll2006l) . GX 1+4 is 
the first and the prototype of the small but growing subclass of 
accreting X-ray pulsars called symbiotic X-ray binaries (SyXB), 
by analogy with symbiotic stars, in which a white dwarf accretes 
from the wind of an M-type gia nt companion (e.g jMasetti et al.l 
I2006L 120071: [Corbet et al.ll2008l) . 

ICutler et al.l ([1986!) proposed an orbital period of about 304 days 
based on variations in the pulse period of the neutron star of 
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GX 1 +4 during the spin-up phase in the 1 970s. Other authors 
(e.g.. iPereira et al.lH999t iBraea et al.ll2Q00h have supported this 
using observations from 1991 to 1998, when the source was al- 
ready in its long-term spin-down phase (see also Fig. [6]). More 
recently, infrared observations have shown a 1161 -day period 
single-line spectroscopic binary orbit, which exclud es the 304- 
day period as the orbital period dHinkle et al.ll2006h . GX 1+4's 
X-ray light curves show strong variability on a timescale of days 
to years but with no modulation on either the optical 1161-day 
orbital period or the previously re ported 304-day X-ray (e.g., 
iNaik et al.ll200l ICorbet et aT1l2008h . 

Based on X-r ay and infrared observations, 

IChakrabartv & Rochd (1 19971) have constrained the distance 
range to the system as between 3 and 15 kpc, depending on the 
evol utionary state of the red giant. For a first-ascent gi ant branch 
star (IChakrabartv & Rochd fl9"97t IHinkle et al.1 120061). the esti- 
mates have narrowed down to 3-6 kpc (IChakrabartv & Rochd 
11997b . The distance estimate by IHinkle et all d2006l) . 4.3 kpc. is 
con sistent with this. Thu s, in spite of its position in the sky (see, 
e.g JPredehl eTail (119951) 1. GX l+4/V2116 0ph is clearly not an 
objec t associated wi t h the c entre of the Milky Way dHinkle et alj 
120061) . IHinkle et al.l d2006l) . however, do not provide any uncer- 
tainty on the distance, and an uncertainty of 0. 1 kpc, based on 
the precision of their estimate, seems to be too precise. They do 
provide an uncertainty in the effective temperature of the giant 
star (±200 K) and its radius (+42,-30 R ), which would give a 
range in the absolute luminosity, and therefore, an uncertainty 
in the distance may be estimated, while taking the interstellar 
extinction into account. Estimates of the latter, in terms of 
E(B - V), have been report ed by vari ous authors: 1.62±0.1 9 
( IChakrabartv & Rochd Il997l). 1 7+0.4 dDavidsen et all 1 19771) 



2.1+0 .1 dJablonski et al.rfl997l) . 2.30+0.06 dShahbaz et all 



|1996|) P1 This rather wide range in extinction has a strong effect 
on the uncertainty in the distance, i.e., about a few kpc. With 
the mean observed V-band magnitude, V= 18.40 + 0.03 mag, 
E(B - V) = 1.62 ± 0.19 o r, equivalently, A v = 5.0 ± 0.6 
IChakrabartv & Rochd dl997l) . and assuming an absolute mag - 
nitude M v of about -0.25 for an M6III star dThe et al.lll990T) . 
we derive a distance of 5.4 + 1.6 kpc. A higher reddening 
value, for exam ple, of E(B - V) — 2 .30 + 0.06 or, equivalently, 
A v = 7.1 + 0.2 dShahbaz etalll 19961) wou ld lead to a distance of 
2 + 0. 2 kpc. We, therefore, conclude that IChakrabartv & Rochd 
d 19971) give a more realistic estimate than IHinkle et alJ d2006l) 
and that the uncertainty on the distance is about 1.5 kpc. In this 
paper we adhere to the distance estimate of 4.3 kpc. 
Currently, there is no well established value for the magnetic 
field of GX 1+4 . Assum i ng the standa rd accretion-disc theory 
dGhosh & Lambl[l979a1lbL IWanl fl987l) . the magn etic field has 
been estimated to be B ~ 10 13 - 10 14 G (e.g., iDotani et alJ 
I1989L iMonv et al.l fT99ll ICui & Smith! l2004h . which is among 
the largest measured for any accreting X-ray pulsar. There have 
been marginal reports of cyclotron scattering resonance rea- 
tures (CRSFs) in the X-ray spectra, which points to a value of 
B ~ 10 12 G for the magnetic field, i.e., u p to two orders lower 
(e.g jRea et aTll2005tlFerrigno et al.ll2007l) . 
GX 1+4's unusual long-term spin behaviour has attracted con- 
siderable interest for many years. Studying the pulse-period evo- 
lution in an accreting X-ray pulsar and relating it to, e.g., the lu- 
minosity changes, allows testing theoretical models and gaining 
insight into the interaction between the pulsar's magnetosphere 



and the accreted matter. In this paper we extend the investiga- 
tion of the spin-period history of GX 1+4 with new observations 
obtained by BeppoSAX, INTEGRAL, and Fermi. The available 
measurements of the pulse period of GX 1+4 span a period of 
about 40 years. This gives us a unique insight in the pulse pe- 
riod evolution of this SyXB. Apart from GX 1+4 still spinning 
down overall, we find irregular trends on top of this evolution. 
We correlate these features with the high-energy flux of the sys- 
tem using the same observations as well as those obtained with 
Swift, and provide an explanation for the spin history seen. 



2. Observations 

2.1. CGRO 

The Compton Gamma-Ray Obse rvatory, CGRO, was a NASA 
mission launched in Ap ril 1991 dGehrels et al.1 [l 9931) and op- 
erative until June 2000 dKaneko et al.l 120061) . The Burst And 
Transient Source Experiment, BATSE, onboard CGRO was a 
sensitive all-sky instrument that consisted of eight uncollimated 
Na I scintillation detectors at the corners of the spacecraft. It cov- 
ered a broad energy range from 15 keV to 100 MeV. Each detec- 
tor module contained a large-area detector (LAD) optimized for 
sensitivity and directional response and a spectroscopy detector 
(SD) optimized for br oad energy coverage and energy resolution 
dFishmanetal.1119921) . 

BATSE monitored pulse frequencies and X-ray pulsed fluxes fo r 
roughly half of the known X-ray pulsars dNelson et al.ll 199Tbh . 
GX 1+4 results from this survey are publicy available in the form 
of pulse source histories covering the energy range from 20 keV 
to 50 keV and daily frequencies from April 1994 to May 19970 



2.2. BeppoSAX 

The X-ray satellite "Satellite per Astronomia X", BeppoSAX, 
was an Italian/Dutch mission launched in April 1996 and oper- 
ated until April 2002, then deorbited in April 2003. It covered 
more than three decades in energy (about 0. 1 to 300 ke V) with 
relatively good energy resolution, and provided imaging capa- 
bilities in the range of 0. 1-10 keV. Together with the Wide-Field 
Camera's (WFCs), the broad-band Narrow-Field Instruments 
(NFIs) provided the opportunity to study the broad-band b e- 
haviour of s everal classes of X -ray sources dBoella et al.ll 19971) . 
The WFCs dJager et al.ll 19971) were two identical coded-aperture 
instruments onboard BeppoSAX. The field of view was 40°x40° 
full width zero response (FWZR), the angular resolution 5' full 
width half maximum (FWHM) and the source-location accu- 
racy was generally better than 1' (99% confidence). The de- 
tectors were sensitive to the energy range 2 to 28 keV. The 
WFCs pointed in opposite directions with respect to each other 
and perpendicular to the NFIs. The pointing directions of the 
WFCs were usually governed by the observations of the NFIs. 
However, a few dedicated campaigns on the Ga lactic bulge re- 
gion were performed by the WFCs (see, e.g., lin't Zand et al.l 
12004 . We analysed WFC data obtained from August 1998 to 
August 2000. 

2.3. INTEGRAL 



1 The observed optical extinction is mainly attributed to interstellar 
dust, while the v ariable X-ray absorption is po ssibly due to the wind of 
the giant star, see IChakrabartv & Rochd JT997h . 



INTEGRAL dWinkler et al.ll2003l) is an ESA scientific mission 
that was launched in October 2002. It is dedicated to spec- 



" See f tp : / /leg acy . gsf c .nasa . gov/ compton/dat a/bat se/| 
pulsar/histories/. 
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Table 1. INTEGRAL observation periods (OP) 



OP 


Start Date 


End Date 


MJD 


Effective 




(UT) 




Exposure 


1 


2005/02/18 


2005/04/19 


53419-53479 


258 ks 


2 


2005/08/16 


2005/10/26 


53598-53669 


377 ks 


3 


2006/08/9 


2006/04/21 


53775-53846 


276 ks 


4 


2006/08/13 


2006/10/23 


53963-54031 


461 ks 


5 


2007/02/15 


2007/04/21 


54146-54211 


853 ks 


6 


2007/08/19 


2007/10/14 


54331-54387 


599 ks 


7 


2008/02/11 


2008/04/20 


54507-54576 


225 ks 



Large Area Telescope (LAT) and the Gamma-ray Burst Monitor 
(GBM). Fermi is dedicated to measuring the cosmic gamma-ray 
flux in the energy range 20 MeV to > 300 GeV, with supporting 
measurements for gamma-ray bursts (Ritz et al. 2009). 
Since 2008 August 12 GX 1+4 ha s been continuously monitored 
by the GBM dMeegan et alJ2009h . The GBM is an all-sky instru- 
ment sensitive to X-rays and gamma rays with energies between 
~8 keV and ~40 MeV. Timing analysis is carried out with chan- 
nels 1 and2oftheNaIdetectorCTIMEdata(12-50keV, 0.256s 
time resolution). 

We analysed GBM data from August 2008 to February 2010. 



troscopy and imaging of celestial y-ray sources in the energy 
range between 15 keV and 10 MeV with simultaneous monitor- 
ing in the X-ray and optical energy ranges. Several instruments 
are on board: SPectrometer onboard INTEGRAL (SPI), Imager 
on Board INTEGRAL Satellite (IBIS), Joint European Monitor 
X-rays (JEM-X) and Optical Monitoring Camera (OMC). The 
former three instruments all collect photons through wide-field 
coded masks. 

IBIS (lUbertini et al.ll2003b comprises two detector planes. The 
field of view is 29°x29° (FWZR) and the angular resolution 12' 
(FWHM). In this paper we only used data collected with one of 
its detectors: the INTEGRAL Soft Gamma-Ray Im ager (ISGRI) 
It is sensitive in the ^15keV to lMeV range dLebrun et alj 
120031) . 

The INTEGRAL Galactic bulge monitoring programme started 
in February 2005 and was initiated to monitor the Galactic bulge 
region on a regular basis in mainly the hard X-ray band. One 
complete hexagonal dither pattern (7 pointings of ^ 1 800 s each) 
is performed during each INTEGRAL orbit around the Earth, 
(roughly every 3 days), whenever the bulge region is visible 
by INTEG RAL: twice per year f or a total period of about four 
months (see lKuulkers et al1l2007l) . 

We analysed the ISGRI data from observations of the Galactic 
bulge monitoring programme between February 2005 and April 
2008. These observations span several periods, see TableQ] Each 
time span is referred to as observation period, or OP for short. 

2.4. Swift 

Swift is a NASA missio n that was launched in November 
2004 dGehrels et alj |2004|) . The Burst Alert Telescope (BAT; 
iBarfhelmv et al.l2005h onboard Swift is a coded-aperture imager 
with a very wide field of view of about 2 steradians, which oper- 
ates in the 15-150keV band. The BAT angular resolution is 22' 
(FWHM). 

The BAT continually monitors the sky with more than about 70% 
of the sky observed on a daily basis. Results from this survey are 
publicly available in the form of light curves covering the 15- 
50keV energy band on two timescales: a single Swift pointing 
(^20 min) and the weighted average for each day0 We used the 
daily average light curve for GX 1+4 to study its long-term hard 
X-ray flux behaviour. 



2.5. Fermi 

The Fermi Gamma-ray Space Telescope, formerly Gamma-ray 
Large Area Space Telescope (GLAST), is a NASA mission that 
was launched in June 2008. It has two instruments on board: the 



3 |http : / /swift . gsf c ■ nasa ■ gov/ docs/ swif t/ result s/| 
transients/ index . html. 



3. Data analysis 

For INTEGRAL/ISGRI, the data were reduced using version 
7 of the INTEGRAL off-line analysis software (OSA 7), 
distributed by the INT EGRAL Science Data Centre (ISDC; 
ICourvoisier et al.ll2003l) . Each data set in a revolution (about 13 
ks long) has been analysed in the energy range between 20 keV 
and 40 keV, following the steps described in the IBIS Analysis 
User Manual (version 6.0). We used the HJight tool to obtain 
light curves with a time bin of 10 s and barycentric correc- 
tion was applied. The source BeppoSAX/WFC flux was recon- 
structed in the 2-25 ke V bandpass with a time resolution of 2 s. 
Period determinations for both the INTEGRAL/ISGRI and 
BeppoSAX/WFC data have been done for segments o f the light 
curve using an epoch-folding analysis dLarssonl ll996). The data 
for such a segment is folded at a number of different test peri- 
ods, and for each period, the x 1 over the resulting pulse profile 
is computed. A best-fit period is determined by fitting a tem- 
plate function describing how x 2 should vary with a test period. 
The template function takes the time sampling of the data into 
account and, in an iterative procedure, the pulse shape of the 
oscillation. Uncertainties are estimated by Monte Carlo simula- 
tions. A set of synthetic pulse light curves with the same time 
sampling and noise level as the data are created and analysed. 
The distribution of determined periods for these simulations is 
used as a measure of period uncertainty. Pulse period determina- 
tion for INTEGRAL/ISGRI data was possible for most hexag- 
onal pattern observations. From a total of 128 hexagonal dither 
patterns it was possible to obtain 121 light curves (in 7 of these 
128 observations the source was below the detection limit of the 
instrument) and to derive 93 pulse periods, which is ~17% of 
the detections. For 28 detections it was not possible to derive the 
pulse period because the source detection significance was lower 
than 14 and no clear maximum in the^ 2 distribution was found 
owing to the noise. 

The analysis of the Fermi/GBM data is complicated by Fermi's 
continuously changing orientation. All intervals of CTIME data 
from the 12 Nal detectors are selected for analysis where the 
high voltage is on, excluding those containing high-voltage tran- 
sients, phosphorescence events, rapid spacecraft slews, South 
Atlantic Anomaly induced transients, electron precipitation 
events, and gamma-ray bursts. Source pulses are then separated 
from the background by fitting the rates in all detectors with 
a background model and subtracting the best-fit model. This 
model includes bright sources and their changing detector re- 
sponses (including Earth occultation steps), along with quadratic 
spline functions that account for the remaining long-term back- 
ground trends. The spline models have statistical constraints on 
the changes in second derivative between spline segments to 
control the model stiffness. These fits are made jointly across 
detectors (with common bright source fluxes) but separately 
for each channel of the CTIME data. The residuals are then 
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Fig.l. Top: GX 1+4 CGRO/BATSE daily light curve in the 
energy range 20-40 keV. Middle: Pulse periods derived from 
CGRO/BATSE data, with the period increasing from top to bot- 
tom. Bottom: Residuals of the periods from a linear fit. 



Fig. 3. Top: GX 1+4 INTEGRAL/ISGRI average flux per point- 
ing in the energy range 18 to 40 keV. Middle: Pulse periods 
derived from INTEGRAL/ISGRI data. Note that the period in- 
creases from top to bottom. Bottom: Residuals from the pulse 
periods using a linear fit. 
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Fig. 2. Top: GX 1+4 CGRO/BATSE daily light curve in the 
energy range 20-40 keV. Middle: Pulse periods derived from 
BeppoSAX/WFC data. Note that the period increases from top 
to bottom. Bottom: Residuals of the periods from a linear fit. 
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Fig. 4. Top: GX 1+4 Swift/BAT daily averaged light curve in 
the energy range 15-50 keV. Middle: Pulse periods derived from 
Fermi/GBM data, with the period increasing from top to bottom. 
Bottom: Residuals from the pulse periods using a linear fit. 



summed over detectors with time-dependent weights that are 
proportional to the predicted (phase-averaged) count rates from 
the source. Short intervals (900s) of these combined residuals are 
then fitted with a constant plus a Fourier expansion to determine 
a pulse profile. The profiles are divided into four-day intervals, 
and the pulse frequency and mean profile are then determined in 
each interval with a search of pulse frequency f or the maximum 
of the Y„ (n = 2) statistic dFinger et al.lll999h . Y„ was formu- 
lated to find a pulse frequency from a series of pulse profiles, 
each represented by a finite Fourier expansion, so it accounts for 
possible frequency-dependent non-Poison noise. 

We did not apply Doppler corrections to the data, since the pro- 
jected semi-major axis has not been measured with X-ray ob- 
servations. The latter is mainly due to the high level of torque 
variability at low frequencies. 



4. Pulse period evolution 

4.1. CGRO 

The CGRO/BATSE data can be found in Fig. [TJ The period in- 
creases in time with an almost linear trend, but with some fluc- 
tuations. A linear fit to these data resulted in a slope of (5.61 ± 
0.03) x 10~ 8 s/s with a Pearson r 2 value of 0.994. The deviations 
from this linear fit are shown in the bottom panel of Fig. Q] In 
this figure it is possible to observe a deviation from the linear 
fit of the spin period evolution around MJD 49700, where the 
neutron star spin-down rate decreases in coincidence with an in- 
crease in the X-ray flux. Overall, the deviations in the linear fit 
shown in the bottom panel of this figure are visually comparable 
to the variations in the X-ray flux shown in the upper panel of 
the same figure. 
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Fig. 5. The discrete correlation function (DCF), i.e., the correla- 
tion coefficient between the 15-50 keV X-ray flux and the pulse 
period change v, as function of the time lag. The minimum near 
zero lag implies an anti-correlation between the pulse frequency 
derivative and X-ray flux. 

4.2. BeppoSAX 

The BeppoSAX/WFC pulse periods are shown in Fig. [2] The 
pulse period continues to increase with time. A linear fit to the 
pulse period as a function of time results in a slope of (1.113 + 
0.017) x 10~ 7 s/s with a Pearson r 2 value of 0.998. The residuals 
from this linear fit are shown in the bottom panel of Fig. [2] The 
CGRO/BATSE X-ray pulsed flux contemporaneous to the WFC 
data is shown in the top panel of this figure. The sparsity of the 
pulse period measurements means we cannot see any clear cor- 
relations between the X-ray flux and the spin-period evolution in 
this data set. 

4.3. INTEGRAL 

The results using the INTEGRAL/ISGRI data are shown in 
Fig. [3] The pulse period increases with time with an almost lin- 
ear trend. We find a slope of (1.031 + 0.007) x 10~ 7 s/s with a 
Pearson r 2 value of 0.995. As can be seen in the bottom panel of 
Fig. [3] the last point at MJD 54541 is only a marginal detection 
of the pulse period. Any linear fit not taking this last data point 
into account does not lead to a significantly different slope. The 
deviations from the linear fit not taking it into account the last 
data point are shown in the bottom of Fig. [3] In this data set 
there is a maximum for the deviations of the spin-period linear 
fit around MJD 53900, followed by a maximum for the X-ray 
flux around MJD 54100. 

4.4. Fermi 

The pulse period results using the Fermi/GBM data can be 
seen in the middle panel of Fig. [4] Again, the period evolu- 
tion for GX 1+4 follows an almost linear trend with a slope 
of (1.0697±0.0002) x 10~ 7 s/s and a Pearson coefficient of 
r 2 = 0.999. The evolution of the Swift/BAT hard X-ray flux is 
shown at the top of this figure. 

4.5. X-ray flux versus spin period 

Theoretical models predict certain correlations between X-ray 
fluxes and pulse period evolution. F or example, the complete 
and continuous da t a series in Fig. [Tlled lChakrabartv et alj d 1997b 
and iNelson et al.l dl997al) to demonstrate for the first time a 



negative correlation between spin-up rate v and X-ray lumi- 
nosity. Motivated by this, we searched for correlations in the 
data shown in Fig. |4] which represents the most complete and 
continuous data series of period and fluxes. We calculated the 
discrete correlation f unction, DCF (see lEdelson & Krolikll 1988*1 
iPeterson et alii 1998b . between v, estimated over eight day bins, 
and the one-day binned Swift/BAT 15-50 keV flux. The DCF for 
the full MJD 54700-55240 time range (Fig. E) shows a strong 
negative correlation at zero lag (-1.3 ±2.6 days). Since the main 
uncertainty of the DCF is random correlation with erratic X-ray 
flickering in the light curve, we divided the data into four sec- 
tions and computed the DCF for each. Strong anticorrelation 
(~ -0.6) was seen in the last two of these, and it was weaker 
(~ -0.3 to -0.4) in the first two sections. In all four the DCF 
minimum for lags between -40 and +40 days occurred close to, 
and was consistent with, no lag. 

Comparing the instantaneous spin frequency derivative derived 
from the Fermi/GBM data with the flux simultaneously mea- 
sured by the Swift/BAT monitor we find a dependence -v oc 
,p(0.30±o.07). This is in line with the correlation of instanta- 
neous spin-down torque with X-ray flux discovered by BATSE 
dChakrabartv et al.lll997h . 

4.6. Long-term pulse period evolution 

The long-term pulse period evolution of GX 1+4 is shown in 
Fig. [6] We combined all measurements that we are aware of 
from the literature with the pulse periods determined from 
BeppoSAX/WFC, INTEGRAL/ISGRI, and Fermi/GBM data as 
presented in the previous subsection. The main features in Fig. [6] 
are the well known switch from a strong spin-up to a spin-down 
trend in the 1980s and a continued spin-down slowly increas- 
ing in average P over the years. The spin-down rate observed 
around 2004 was P ~ 10~ 8 s/s, while the spin-down trend ob- 
served in this work with INTEGRAL/ISGRI (2005-2008) and 
with Fermi/GBM (2009-2010) is P ~ 10~ 7 s/s, taking the spin- 
down measure in frequencies, the global evolution since the 
start of spin-down is described very well by a linear trend of 
-0. 1 177(3) mHZ/y as it is shown in the upper panel of Fig. Q 
On top of these spin-down trends, irregularities are seen that 
have sometimes been proposed to correlate with the binary orbit 
(e.g.. |Pereira et al .ll999HBraga et al.l2000h . (For a more detailed 
description see Section [57X1 and Appendix lAl). These irregulari- 
ties are reflected in the middle panel of Fig. [7] but we note also 
that other apparent trend changes occur at times far from the pre- 
dicted perigee passages. 

5. Discussion 

5.1. A retrograde disc in GX 1+4? 

When GX 1+4 was discovered, the pulsar was spinning up, and 
the spin-up rate v was apparently correlated with the observed 
X-ray flux, which is usually interpreted as a direct measure of 
X-ray luminosity, hence m ass accretion rate M (e.g jDotv et al.l 
ll98l[lRicketts et alJ[l982l) . This supported the idea that the pe- 
riod decrease is produced by ac cretion torques of a prograde 
disc a round the neutron star (e.g jGhosh & Lamblll979bl:IWand 
Il987h . 

Arou nd 1983-1984 GX 1+4 entered a lo w X-ray luminosity 
state (lHall & Davelaari[T98l lMukail[l988h . suggesting a large 
reduction in the mass ac cretion rate, and it started to spin down. 
iMakishima et al.l (11988) were the first to propose for GX 1+4 
that a retrograde disc is formed around the neutron star by matter 
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Fig. 6. The long-term pulse period evolution of GX 1+4, incorporating data from the literature and results of this work. See Table lB.il 
for data and references. 



captured from the stellar wind of the M-giant companion of the 
pulsar, to explain the observed spin-down. Since then, this ret- 
rogra de disc hypothesis has been supported by several authors 
(e.g.. IDotani et al.|[l989t IChakrabartv et al.|[l997t iNelson et ail 
Il997al) for different reasons. One of the arguments in favour of 
the retrograde disc hypothesis is related to the magnetic field. 
To explain the spin-down of a pulsar accreting from a pr o- 
grade stand ard disc (see e.g., IGhosh & Lambl I1977L Il979allbl 
IWangl Il987h . the neutron star must be rotating very near its 
equilibrium spin period. The pulsar should then be in a quasi- 
equilibrium state, and GX 1+4 consequently has the strongest 
kno wn magnetic field of any neutron star (B ~ 10 13 - 10 14 G 
(e.gjDotani et alJl989l:lMonv et al.ll99"TtlGreenhill et al.ll993t 
ICui& Smith] 120041) . The retrograde disc scenario would elim- 
inate the need for such an unusually strong magnetic field, 
but the pulsar wo uld necesarily be far f r om its equilibrium 
spin period (e.g.. [Makishima et all [19881: IDotani et aD Il989t 
IChakrabartv et al.lll997UNelson et al.lll997al) . 

The negative correlation of spin-up rat e v and X-ray flux found 
in BATSE obser vations in the 1990s (IChakrabartv et alj|1997t 
iPaul et al.lfT997b and in this work is also not consistent with 
the standard disc accretion model by IGhosh & Lambl (Il979bl) . 
which predicts higher spin-up rates for higher X-ray luminosi- 
ties. However, it is consistent with a retrograde disc around the 
neutron star, taking angular momentum off the pulsar, so that 
higher spin-down rates for higher X-ray luminosities are ex- 
pected. However, it remains a chief question how a retrograde 



disc could form and remain stable over such a long time period 
(about 30 years). 

On the other hand, positive correlations between spin-up rate v 
and X-ray luminosity have also been found du ring the steady 
spin-down trend (e.g.. IChakrabartv et al.l [1991 . In addition to 
these positive correlations, during this spin-down steady trend 
some spin-up episodes related to bright flares, where the source 
has reached almost the same luminosity as it had in the 1970s, 
have taken place (e.g., IChakrabartv et alj|1997t iFerrigno et alj 
l2007h . Again, it is not possible to explain this with accretion 
from a retrograde disc. 

The presence of fast flickering has been interpreted as evi- 
dence of acc retion discs around neutron stars in general (e.g., 
lHorn3[l99l and in particular, for GX 1+4 on t he basis of 
such flickering in the opti cal light curves (e.g., iBraga et all 
[1991: Ijabionski et alJ [1991 . Meanwhile, the absence of this 
flickering when the system is faint h as been interpreted a s 
episodic interruptions of the accretion dJablonski et al.l [T99l . 
Fu rthermore, coherent op tical pulsations of two minutes found 
bv lJablonski et al.l d!99'7 



in an accretion disc ( 


Chested 


Chakrabarty & Roche 


119971). 



Moreover. IChakrabartv & Roche! ( 1 19971) have carried out an op- 
tical emission-line diagnostic study of the optical and infrared 
spectra, which also suggested there is an accretion disc around 
the neutron star. 

According to numerical simulations of mass accretion onto 
the neutron star in wind-fed systems, the accreted specific 
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by BATSE, RXTE/ASM, and Swift/BAT, aver- 
aged over 30-day intervals. 



angular momentum could change sign in an erratic manner, 
whic h may lead to alternating spin-up and spin-down episodes 
(e.g. iTaam & Frvxell|[T988t iMatsuda et all 19911: iFoglizzo et all 
120051) . This phenomena, known as " flip-flop ", has been used 
to explain the random walk of the spin period of sources like 
Vela X-l (e.g jBovnton et aDll984l) . However, it is important to 
note that these sim ulations apply to supergi ant stars with power- 
ful winds (see e.g. iKudritzki & Pu is 2000), and not to M-giant 
starsQ 

Last but not least, the X-ray luminosity increase accompanied by 
a rather constant sp in-down rate observed by, e.g., Ginga in 1987 
dSakao et al.lll990l) implies (temporary) luminosity-independent 
spin-down behaviour. This cannot be understood in terms of ei- 
ther the prograde disc, or the retrograde one. 

To conclude, ever since GX 1+4 entered the low-luminosity state 
around 1984, overall it kept spinning down with the general spin- 
down accompanied by spin-up episodes during flares. The pres- 
ence of both positive and negative correlations between X-ray 
flux and spin-down rate, and episodes of no correlations between 
X-ray flux and spin-down rate, are not possible to explain either 
with standard disc accretion or with a retrograde disc. One can- 
not explain the long-term steady spin-down with an alternating 
disc either, because the neutron star tends to spin-down (except 
for short episodes). It can also not be explained by a retrograde 



disc, as the spin-down rate is increasing with time while the aver- 
age X-ray luminosity0 is stable around Lx ~ 10 35 - 10 36 erg s _1 . 



5.2. Quasi-spherical accretion in GX 1+4 

In this section we discuss how the long-term spin-down be- 
haviour of GX 1+4 presented in Fig.|6]can be explained by quasi- 
spherical accretion onto the neutron star from the stellar wind 
from the secondary companion. 

In wind-fed pulsars with long orbital periods, accretion onto the 
neutron star can proceed quasi-spherically; i.e. an accretion disc 
around the neutron star magnetosphere cannot be formed at all. 
The formation of the disc depends on whether the specific angu- 
lar momentum of matter j,„ near the magnetospheric radius R,„ 
is larger or smaller than the Keplerian value jic(Rm) - y/GMR m . 
Assuming the specific angular momentum conservation, j m can 
be related to that of gravitationally captured stellar wind matter 
in the zone of bow shock at the Bondi radius Rq = 2GM/(Vf v + 
V 2 rb ) ~ 10 12 v^ 2 cm (here v-/ is the relative stellar wind velocity 
in units of 100 km/s) , j,„ = j w . To within a numerical factor 



of order one, 



./>.■ 



(27T/P orh )Rl * 3 x 10 I7 (iW1000'') 



d\-l„-4 



which can be smaller than jxiRm) for typical magnetospheric 
radiS R m = 10 cm, so the formation of quasi-spherical accre- 
tion flow is very likely. 



4 V2116 Oph is an M6-giant l lChakrabartv & Rochd 119971 : 
IHinkle et all l2006h . Such stars ha ve slow and strongly variable 
stellar winds (see e.g., Crowlev 2006). 



5 The luminosity value has been corrected to the distance of 4.3 kpc 
l lHinkle et all20061) . as are the luminosity values used hereafter. 

6 The value of the magnetospheric radius is determined by the pres- 
sure balance at the magnetospheric boundary, and generally can be no- 
tably different for disc and spherical accretion 
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There can be two very different regimes. If the X-ray lumi- 
nosity of the central source is high enough (~ 10 37 erg/s), the 
Compton cooling of plasma in the region of the bow shock 
is rapid, so the matter freely (supersonically) falls toward the 
magnetosphere, near which the shock is formed. This regime 
with different deg rees of physical desc ri ption was con s idered , 



1983), 



1990). 



e.g., in papers by lArons & Leal ( 1976h. Burnard et alJ 
iBisnovatvi-Kogan] d 19911) . or llllarionov & Kompaneetsl 
In contrast, the cooling can be ineffective for moderate and low 
X-ray luminosities, and the matter moves toward the neutron 
star magnetosphere subsonically, by forming a hot quasi-static 
shell around the magnet osphere (the settling acc retion). Such 
shells were considered bv lDavies & Pringle1(ll980h . A superadi- 
abatic temperature gradient can be established in the shell, lead- 
ing to development of large-scale convective motions and turbu- 
le nt cascades. Recently, thi s regime of accretion ha s been studied 
by Shakura et alJ din press!) (see also first results in lPostnov et alJ 
120101) . In the settling accretion regime, the accretion rate onto 
the neutron star is determined by the ability of plasma to en- 
ter the magnetosphere via instabilities, and this determines the 
mean radial velocity of matter in the shell by the mass conser- 
vation law. In the last paper it was found that the critical X-ray 
luminosity below which the settling accretion regime sets in is 
about 4 x 10 36 erg/s. 

In the free-fall accretion regime, the neutron star spin behaviour 
is determined by the sign of the specific angular momentum 
of captured matter (prograde or retrograde). Nevertheless, in 
the case of GX 1+4, this regime seems to be unlikely for two 
reasons: (a) the X-ray luminosity of the source at the spin- 
down stage is quite low and (b) the long-term spin-down with 
short episodic spin-ups are difficult to reconcile with expected 
features of wind accretion - although alternation of the pro- 
grade/retrograde angular momentum of the captured matter in 
inhomogeneous stellar wind is possible, the predominance of the 
retrograde sign is very enigmatic. 

In contrast, in the settling accretion regime, the hot convective 
shell can mediate the angular momentum transfer to/from the 
neutron star magnetosphere, and the neutron star can spin up or 
down depending on the sign of the difference of angular veloc- 
ity of matter near the magnetospheric boundary and that of the 
magnetosphere itself. The gas-dynamic problem of spherically- 
symmetric accretion flow with cooling and heating due to tur- 
bulence (generally, anisotropic) was considered in Shakura et al. 
(2011). It was found that in this regime with increasing X-ray 
luminosity, the neutron star spin-down can change for spin-up 
(even abruptly, if L x exceeds the critical value for the existence 
of the shell) and vice versa, and the fluctuations in the spin fre- 
quency can anti-correlate with flux fluctuations, which indeed 
follows from the analysis of BATSE and Fermi/GBM observa- 
tions of GX 1+4 (see above). 

It therefore seems very likely that the settling accretion is un- 
der way at the present low-luminosity spin-down state of GX 
1+4. Then the source must have had higher X-ray luminosities 
in the preceding long-term spin-up state. Indeed, at high X-ray 
luminosities, the Compton cooling near the magnetosphere is 
very strong and a free-fall gap in the matter flow appears above 
the magnetosphere, the casual connection between the magneto- 
sphere and the shell is lost, and only spin-up of neutron star is 
possible. 

The fluxes reported during the long-term spin-up phase prior to 
the 1980s and the luminosities derived from these are generally 
a few times higher than those reported during the still ongoing 
spin-down phase. Accounting for the differences in instrumenta- 
tion and energy ranges and the corresponding uncertainties, we 



arrive at l uminosities L2-60 kev ~ 1 - 5 x 10 37 erg s _1 for the early 
data (e.g jDotv et al.llT~98 lh . After the torque reversal, luminosi- 
ties in this energy interval mostly remain in the range of a few 
times 10 36 erg s _1 . A flare seen by BATSE dChakrabartv et al.l 
119971) had a pulsed flux corresponding to L2-6O kev ~ 9 x 10 
erg s _1 , which indicates a comparable brightness to the early 
data. We refer to AppendixlAlfor a more detailed discussion. 
The middle panel of Fig. [7j shows long-term quasi-periodic 
frequency fluctuations. Their fractional amplitude A<w/|<jj s d| = 
(6co AT) =i +0.1 over a time interval of about one orbital period 
(At 1.5 x 10 8 s) corresponds to the fractional derivative ra- 
tio to/|di s dl — 0.3. Apparently, they are marginally correlated 
with periastron passages of the binary system, but not clearly 
correlated with the X-ray flux variations indicated in the bot- 
tom panel of Fig. Q These long-term fluctuations can be due to 
smooth variations of the wind density p w and velocity near the 
gravitational capture radius. 

The short-term spin-up episodes sometimes observed on top of 
a stea dy spin-down behaviour (see Fig. 2 in IChakrabartv et al.l 
1997, and Fig. [TJof this paper) are correlated with an enhance- 
ment of the X-ray flux, in contrast to the negative frequency-flux 
correlation at the spin-down discussed above (see Figs.|4]and|5]l. 
During these short spin-ups, at was about half of the average 
d» su observed during the steady spin-up state of GX 1+4. The 
X-ray luminosity during these episodic spin-ups was approxi- 
mately five times higher than the mean X-ray luminosity dur- 
ing the steady spin-down. These facts are consistent with the 
quasi-spherical accretion model that predicts transitions from 
spin-down to spin-up with increasing mean accretion rate and 
the reversal of flux-period correlation properties (see Shakura et 
al. 2011 for more detail). However, an increase in the accretion 
rate by more than one order of magnitude in GX 1+4 could de- 
stroy the shell due to rapid radiation cooling and the establishing 
of the free-fall accretion regime onto the magnetosphere. In that 
case only spin-up of neutron star is possible. 

6. Conclusions 

We present the results of observing the pulse period behaviour 
of the symbiotic X-ray binary GX 1+4. New measurements 
by BeppoSAX, INTEGRAL, and Fermi confirm the continuing 
overall spin-down (about ~ 3.4 s/y) of the neutron star since the 
torque reversal in the early 1980s. The pulse period has increased 
by about ~5Q% over the past three decades and has reached the 
highest value ever observed for this source. The X-ray luminos- 
ity during the extended spin-down phase has in general been sig- 
nificantly lower than during the spin-up phase of the 1970s (see 
AppendixlAl. 

The global spin-down follows a linear trend in angular frequency 
with deviations < 2% (see Fig. [7]). During this time interval, 
there ha ve only been brief insta nces of spin-up observed by 
BATSE dChakrabartv et alj|1997l) related to bright flares where 
the X-ray luminosity has almost reached the value it had in the 
1970s, and a possible spin-up in 2004 observed by INTEGRAL 
dFerrigno et al.ll2007l) . 

On top of the long-term spin-up and spin-down trends going on 
for decades or more independently of the X-ray luminosity, there 
are pulse period fluctuations on shorter timescales that do show 
anti-correlation with the source X-ray flux as demonstrated in 
Section|4]and, apparently, marginally correlated with the orbital 
phase. 

As shown in Secs l5.ll and l5.21 these observational facts are dif- 
ficult to reconcile with prograde or retrograde disc accretion but 
can be explained by assuming quasi-spherical accretion onto the 
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neutron star from the stellar wind of the M-type giant compan- 
ion. 
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Appendix A: Historical spin period and X-ray flux evolution of GX 1+4 

In this section we attempt to summarize the evolution of the spin period, as well as of the brightness of GX 1+4. For a table of all 
pulse period measurements, we refer to Appendix B. The former is rather straightforward, the latter much less so. The difficulty 
arises for various reasons. First of all, the flux values come from a variety of instruments with different energy ranges and spectral 
responses. Second, few public ations give enoug h detail to determine the exact spectral shape used in the analysis. Finally, different 

broad-band observations (e.g. jNaik et al.|[2005h do demonstrate clear spectral variatio n, including a varia ble absorption. 

To get a basis for comparison of flux we used the detailed spectral shapes published by Naik et af] (120051) and lFerrigno et al.l (120071) 
to simulate a range of spectra under different conditions. We varied crucial spectral parameters within the uncertainties given for 
these spectra to cover more of the possible parameter range. From these simulated spectra we then derived simulated fluxes in 
various energy ranges and used the ratio of these fluxes to derive flux conversion factors between the bands, together with the 
correspondi ng uncertainties. Ba sed o n this analysis , we co nverted selected fluxes to the 2-60 keV band, for which fluxes were 
reported b y ICui & Smith! (|2004|) and iFerrigno et al.l d2007l) and very close to the SAS 3 band (1.2-55 keV) for which fluxes are 



reported in lDotv et al.l d 198 lh . 



The uncertainty in such conversions can be fairly high, especially when going from a narrow band to a broader band. In our analysis 
we estimate these to be up to, e.g, ~55% scaling a 2-6 keV flux to the 2-60 keV band, while allowing for uncertainties in the soft 
band spectral shape and the level of absorption at any given time. Therefore, we concentrate as much as possible on fluxes that can 

be compared directly for our discussion. 

For an earlier overview of the flux evolution of GX 1+4 between 1970 and 1988 we refer to Table 2 of iMcClintock & Levenfhail 
(fl989l 



A. 1. Early observation - the spin-up phase 

After its detection in the 1970s GX 1 +4 was spinning up with the fastest rate (P 8 .5 x 10~ 8 ss~' or -2.55 syr~') among the 

known X-ray pulsars at the time (e.g. jDotv et al.l[T98lHElsner et al ] U985HNagaselll989h . The OSO-7 scans between October 1971 
and May 1973 found flux values in the Argon counters varying between ~ 2 and ~ 12 cts s _1 dMarkert et al.ll 1979b . whi ch roughly 
scale J] to a flux in the 2-60 keV band of (3 + 1 - 18 + 6) x 10~ 9 ergs cirT 2 s" 1 . From SAS3 observations in 1975 and 1976. lDotv et al.l 
(1 19811) found fluxes i n the 1 .2-55 keV ra nge of (4.5 +1 - 22 + 10) x 10~ 9 ergs cirT 2 s _1 , comparable, within the uncertainties, to 
OSO-8 observations (iBecker et al.lll976l) . 



A.2. Low state and torque reversal 

EXOSAT failed to dete ct the source in 1983 a nd 1984, revealing an extend ed low state w ith upper limits on the 2-10 keV flux of 
< 10 ergs cm 2 s 1 dHail & Davelaar Il983h and < 10 11 ergs cm 2 s 1 dMukaill 19881) . respectively. Even allowing for a factor 
3-7 when scaling this to the broader energy range above and accounting for various uncertainties, the source was fainter by two 
orders of magnitude or more during the spin-u p phase. 

A first detection by Ginga in March 1987 dMakishima et alJ 119881: iDotani et all Il989l) found the source still very faint at 
~ 10~ 10 ergs crrT 2 s _1 in the 2-20 keV band (sca ling to roughly twice this value in 2-60 keV). The pulse period of 1 10.223 + 0.003 s 
was longer than the last previous measurement dRicketts et al.ll 1982b . demonstrating that a torque reversal had taken place and the 
source was spinning down (see Fig. [6] and Table IB. il l. In addition, the pulse profile of that early observation was peculiar and very 
different from previously observed ones. 



A. 3. Early spin-down phase 

HEXE and Ginga observations between 1987 and 1989 (IDotani etal.ll 19891: ISakao et all 19901: iMonv et al.ll 19911) revealed a remark- 
ably constant spin-down rate at P ~ 4.5 x 10~ 8 s s in spite of a further overall increase in the X-ray bright ness. Ginga fluxes from 
March 1988 and August 1989 were reported at five to six times that of the first detection in March 1987 dSakao et al.ll 1990h . No 
correlation between spin torque and X-ray flux was reported during this period. 

GX 1+4 was observed by GRANAT/SIGMA during 1990 and 1991 (iDenis et al.lll991b . The spin-down of the pulsar was confirmed, 
but at a lower rate (P ~ 2.6x 10~ 8 s s~'). The linear spin-down trend was recovered around September 1991 with a slightly increased 
rate over the one observed before 1990 dMandrou et al.l[l994l) . 

Ginga observations in September 1990 and September 1991 showed a decrease in the 2-20 keV X-r ay fluxes to ~ 6 a nd ~ 2 x 
10~ 10 ergs cirT 2 s _1 , respectively, suggesting another low state, together with a drastic increase in N R dKotani et al.lll999h . Despite 
this, the spin-down rate was found to be very stable at P ~ (3 - 4) x 10~ 8 s s _1 . 



A. 4. Regular monitoring - spin-down with interruptions 

BATSE onboard CGRO monitored GX 1+4 daily between April 1991 and September 1994. The average X- ray flux in the 20-60 keV 
range was ~ 2 X 10~ 10 ergs cirT 2 s _1 , but interrupted by intermittent bright flares of about 20 day duration (IChakrabartv et al.ll 1997b . 
Duri ng this time a negative correlation between spin-up torque and X -ray flux, i.e., stronger spin-down for higher flux, was observed 
(e.g JPaul et all 1997t IChakrabartv et alJll997tlNelson et alj|1997al) . 



7 See also http : / / www . astronomy cafe .net/qadir/ql723 . html 
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In September 1994 an ASCA observation found the source brighter again with a 2-20 keV flux of ~ 9 x 10~ 10 ergs cirT 2 s _1 
dKotani et al.l [l999'). One month later, BATS E detected a ~200 day b right state, when the pulsar started to spin-up, resuming the 
spin-down when the low state was recovered dChakrabartv et al.ll 19971) . 

BATSE observations between May 1995 and June 1996 showed an erratically varying X-ray luminosity, although the pulsar contin- 
ued to spin-down at a rel atively steady rate with a weak negative peak at zero lag in the cross-correlation of torque and flux histories 
dChakrabartv et al.l 1 997|) . After a period of low activity for about a month, GX 1+4 reappeared briefly in a bright flare, during which 
it transited to spin-up ( Chakrabartv et al.ll 19971) . The pulsed flux in this flare was ~ 2 x 10~ 9 ergs cm 2 s _1 in the 20-60 keV band. 
According to our analysis, this would correspond to a ~ (3 — 4.5) x 10~ 9 ergs cirT 2 s _1 pulsed flux in the 2-60 keV broad band. 
About ten days after this flare, the source dropped below the the detection limit and remained undetected until December 1996. 
During this low state, a series of RXTE observations also failed to detect GX 1+4 , giving an upper limit for the X-ray luminosity 
that was comparabl e to the 1983/1984 EXOSAT low state dChakrabartv et al.lll997l) . 

ICui & Smithl d2004l) monitored GX 1+4 weekly with RXTE in 2001 and 2002. While the source was "quite bright" at the beginning 
of the year 2002 with fluxes (2-60 keV) of ~ (2 - 3) x 10~ 9 ergs cirT 2 s _1 , it made a transition to a faint state in mid June, reaching 
a minimum flux of ~ 3 x 10~ n ergs cirT 2 s _1 , again comparable to the EXOSAT upper limits during the time of the torque reversal. 
At fluxes below ~ 2 x 10~ 10 ergs cirT 2 s _1 some observations had no detectable X-ray pulsations, while others presented a clearly 
pulsed signal. Towards the end of the year, the X-ray luminosity i ncreased again. A joint Chandra & RXTE observation in August 
2002 found a 2-20 keV X-ray flux of 3.8 x 10 ergs cm - 2 s" 1 dPaul et al. 2005}. The average spin-down rate during this period 
was P ~ 4.4 x 10~ 8 s s 1 dCui & Smifhll2004l) . 

During INTEGRAL observations between March 2003 and September 2004 dFerrigno et alJl2007h the source increased its flux 
by about a factor of 5, with erratic variations of about one order of mag nitude. For two time periods with broad-band coverage 
by the INTEGRAL instruments in February /March and September 2004, iFerrigno et al.l d2007l) determined the 2-60 fluxes to be 
~ 1.7 X 10~ 9 ergs cirT 2 s _1 and ~ 2.3 x 10~ 9 ergs cirT 2 s , respectively. 

The secular spin-down trend of GX 1+4 at the time was P ~ 6.6 x 10~ 8 s/s, with a spin-up of the source observed during the 
high-luminosity state in September 2004 dFerrigno et alj[2007h . 
For more recent data we refer to Secs|2]and|4] 
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Appendix B: GX 1+4 pulse period measurements 

For reference we have collected all pulse period measurements presented in Fig.[6]and give them in Table lB.fl below, 
Table B.l. Pulse period measurements of GX 1+4 



MJD 


spin 00 




Instrument 


MJD 






Instrument 


ACiQTZ Z 

4Uo / j.j 


135.0 ±4.0 


Mil balloon 


A 1 S.H 1 ^ 
41 J / l.J 


129.0 ±0.3 


Copernicus 2 


A 1 ZHQ Z 
41 J /O. J 


129.6 ± 0.6 


Copernicus 


A 1 Z 
41 /DO. J 


131.7 ±4.2 


Copernicus 2 


4Z1 jy. j 


128.1+0.3 


r>;„ kill™™ 3 

Kice balloon 


4Z0 /{J.J 


122.46 ± 


0.03 


A< A O 4 

UoU-8 


4Zo>'J.Z 


122.34 ± 


0.06 


oAS-3 


4Z / /Z.o 


121.367 d 


: 0.004 


C A C 1 5 


4Zolz. 1 


120.6589 


± 0.0003 


C A C "3 5 

SAS-3 


4Zo J 3. J 


120.493 d 


: 0.003 


OSO-8 


4ZVUV.J 


120.5 ± 0.5 


MOT U*.ll~~~ 7 

JNKL balloon 


4Z9DJ.Z 


120.19 + 


0.05 




4JZ1 j.j 


118.873 i 


: 0.005 


OaO-8 


A1^\ A A C 

4jZ44.j 


118.715 d 


: 0.005 


Ariel 5 8 


43470.6 


117.45 ± 


0.20 


JNKL balloon 


43593.5 


116.49 + 


0.10 


OSO-8 


4j / /Z.U 


1 14.254 d 


: 0.011 


OSO-8 


43834.5 


113.75 + 


0.25 


JV1PI/ luebingen balloon 


43952.5 


112.68 + 


0.03 


binstein/MPC 


44(JoZ.9j 


1 12.076 d 


: 0.003 


A „" 1 C 11 

Ariel 6 


44346.95 


109.668 d 


: 0.003 


A ■ 1 C 11 

Ariel 6 


/i c r\/~ r\ r\ 

45060.0 


108.8 ± 0.2 


Balloon experiment 12 


46052.0 


108.1+2.6 


rrTrn i n 1^ 14 

TIFR balloon • 


46754.25 


111.5 ± 1.2 


Tasmania balloon 15 


46885.0 


110.233 i 


: 0.003 


t-*: 16 

Ginga 


47099.54 


111.03 ± 


0.02 


» «■ ttr i/iirvr 17 

Mir/Kvant/HEXE 


AT~\ A 1- ! 1 

47247.3 


111.59 + 


0.02 


r~<- 18 

Ginga 


A 1 "! A , ~\r\ f 1 

47429.61 


112.359 d 


: 0.005 


a if' itt~ i/iirvr 17 

Mir/Kvant/HEXE 


47765.47 


113.626 d 


: 0.002 


/-<■ 19 

Ginga 


48126.8 


114.540 d 


: 0.062 


/- ~> TA a XT AT"/ A TIT TA ''O 

GRANAT/ART-P 


48171.7 


114.657 d 


: 0.014 


ta a xt a rp / a i~a t ta 7f) 

GRANAT/ART-P 


48337.87 


115.06 ± 


0.03 


-i TA A XT A nn /OTrH It A ''I 

GRANAT/SIGMA 


48347.52 


115.086 d 


: 0.005 


GRANAT/SIGMA - 


48368.5 


115.2280 


± 0.0016 


ri /■ a ivt at" or 72 

CGRO/BATSE 


48373.5 


115.2617 


± 0.0022 


CGRO/BATSE 


48378.5 


115.2932 


± 0.0025 


/"■'f "' fa a iTt AT^or 22 

CGRO/BATSE 


48383.5 


115.3301 


± 0.0016 


CGRO/BATSE 


/I o o o o c 

48388.5 


115.3664 


± 0.0016 


~> TA /- A ITt AT"Or 22 

CGRO/BATSE 


48393.5 


115.4112 


± 0.0014 


CGRO/BATSE 


48398.5 


115.4567 


± 0.0015 


/"*• f ~* ti f ~\ iTt at" or 72 

CGRO/BATSE 


48403.5 


115.5001 


± 0.0013 


/"^/" , n/'A ivt tTCn 2^ 

CGRO/BATSE 


48408.5 


115.5352 


± 0.0013 


/"■'f "' n a /r^ AT^or 27 

CGRO/BATSE 


48413.5 


115.5668 


± 0.0017 


/"i/^ta*"* /ta ATor 22 

CGRO/BATSE 


4 O A 1 O C 

48418.5 


115.5978 


± 0.0025 


~> TA /- A ITt AT"Or 22 

CGRO/BATSE 


48423.5 


115.6323 


± 0.0021 


CGRO/BATSE 


48428.5 


115.6617 


± 0.0018 


z' ~> TA /- A ITt AT" Or 77 

CGRO/BATSE 


48433.5 


115.7029 


± 0.0021 


/"^/" , n/'A / ri tTdn 2^ 

CGRO/BATSE 


A O A 1 O C 

48438.5 


115.7378 


± 0.0017 


/"■'f "' n a /r^ AT^or 27 

CGRO/BATSE 


48443.5 


115.7689 


± 0.0020 


/"i/" 1 !-**"* /ta ATor 22 

CGRO/BATSE 


/in /I in 

48448.5 


115.8091 


± 0.0023 


~> TA /- A /T"^ AT"Or 22 

CGRO/BATSE 


48453.5 


115.8482 


±0.0017 


L(jKU/dA1 bb 


4o4jo.j 


115.8836 


± 0.0020 


CCjKU/dAI bb 


48463.5 


115.9148 


± 0.0025 


/"^/" , n/'A / ri ATdn 2^ 

CGRO/BATSE 


48468.5 


115.9493 


± 0.0018 


/"■'f "' n a /r^ AT^or 27 

CGRO/BATSE 


48473.5 


115.9780 


± 0.0024 


/"i/" 1 !-**"* /ta ATor 22 

CGRO/BATSE 


48478.5 


116.0134 


± 0.0023 


~> TA /- A /T"^ AT"Or 22 

CGRO/BATSE 


48483.5 


116.0383 


± 0.0061 


L(jKU/dA1 bb 


/I O A O O C 

4o4oo.j 


116.0648 


± 0.0052 


CCjKU/dAI bb 


48493.5 


116.0820 


± 0.0038 


/"^/" , n/'A / ri ATdn 2^ 

CGRO/BATSE 


48498.5 


116.1130 


±0.0031 


/"■'f "' n a /r^ AT^or 27 

CGRO/BATSE 


48503.5 


116.1392 


± 0.0040 


/~i ta /" "\ /ta a ry or 22 

CGRO/BATSE 


A O CCAO C 

48508.5 


116.1592 


± 0.0027 


~> TA /- A /T"^ AT"Or 22 

CGRO/BATSE 


48513.5 


116.1799 


± 0.0022 


L(jKU/dA1 bb 


48533.5 


116.3084 


± 0.0044 


CCjKU/dAI 5b 


48538.5 


116.3213 


± 0.0028 


/~1 TA /" A /TA A T 1 or 22 

CGRO/BATSE 


48543.5 


116.3713 


± 0.0014 


/~> TA A /TA A T" O r 77 

CGRO/BATSE 


A O C /I O C 

48548.5 


116.4208 


± 0.0014 


/- • s ■> n /■ a / ATor 22 

CGRO/BATSE 


48553.5 


116.4638 


± 0.0040 


~> TA /- A /T"^ AT"Or 22 

CGRO/BATSE 


48558.5 


116.4749 


± 0.0052 


L(jKU/dA1 bb 


48563.5 


116.5092 


± 0.0029 


/"■/'~>r)/'^/r) AT'or 27 

CCjKU/dAI bfc 


48568.5 


116.5421 


± 0.0034 


/~i ta f a /ta a t 1 or 22 

CGRO/BATSE 


48573.5 


116.5774 


± 0.0014 


S~ ^ TA A /TA A T" O r 77 

CGRO/BATSE 


48578.5 


116.6152 


± 0.0017 


/-• s ■> n /■ a / AT"or 22 

CGRO/BATSE 


48583.5 


116.6493 


± 0.0019 


/--• r~< n /■ a /r^ AT'or 22 

CGRO/BATSE 


A O C O O C 

48588.5 


116.6926 


± 0.0041 


/"■/'~>r>/'A/r) AT>or 2P 

CCjKU/dAI bb 


48593.5 


116.7240 


± 0.0028 


CCjKU/dAI bb 


48598.5 


116.7618 


± 0.0025 


/~i t~\ i~» /" a /n a t 1 or 22 

CGRO/BATSE 


48603.5 


116.7910 


± 0.0019 


/"> ^ TA A /TA A T" O r 77 

CGRO/BATSE 


48608.5 


116.8300 


± 0.0020 


/'•/'■> n /■ a /r^ AT"or 22 

CGRO/BATSE 


48613.5 


116.8659 


± 0.0020 


~> TA /- A /T"^ AT"Or 22 

CGRO/BATSE 


48618.5 


116.9093 


± 0.0016 


/'•/'-> r~» a /n at* or 22 

CGRO/BATSE 


48623.5 


116.9486 


± 0.0017 


/"■'f "' fa /" a /r^ ATor 27 

CGRO/BATSE 


A Of n O C 

48628.5 


116.9918 


± 0.0014 


/~> i~» f a /n a t 1 or ^2 

CGRO/BATSE 


48633.5 


117.0330 


± 0.0012 


/~> TA A /TA A T" O r 77 

CGRO/BATSE 


48638.5 


117.0752 


± 0.0013 


n f a /"n a T" or ^2 

CGRO/BATSE 


48643.5 


117.1177 


± 0.0014 


~> TA /- A /T"^ ATOr 22 

CGRO/BATSE 


48648.5 


117.1563 


± 0.0017 


/'•/'-> r~» a /n at* or 22 

CGRO/BATSE 


48653.5 


117.2024 


± 0.0019 


/"■'f "' fa /" a /r^ ATor 27 

CGRO/BATSE 


48658.5 


117.2400 


± 0.0020 


/~1 z" 1 1~» /" A ITi A T 1 Or ^2 

CGRO/BATSE 


48663.5 


117.2894 


± 0.0014 


/~> TA A /TA A T" O r 77 

CGRO/BATSE 


48668.5 


117.3327 


±0.0011 


CGRO/BATSE 22 


48671.0 


117.3590 


± 0.0160 


GRANAT/SIGMA 23 


48673.5 


117.3786 


± 0.0013 


CGRO/BATSE 22 


48678.5 


117.4225 


± 0.0014 


CGRO/BATSE 22 


48683.5 


117.4645 


± 0.0015 


CGRO/BATSE 22 


48688.5 


117.4987 


± 0.0030 


CGRO/BATSE 22 


48693.5 


117.5280 


± 0.0053 


CGRO/BATSE 22 


48703.5 


117.6059 


± 0.0027 


CGRO/BATSE 22 


48708.5 


117.6299 


± 0.0041 


CGRO/BATSE 22 


48713.5 


117.6581 


± 0.0036 


CGRO/BATSE 22 


48718.5 


117.6943 


± 0.0024 


CGRO/BATSE 22 


48723.5 


117.7281 


± 0.0034 


CGRO/BATSE 22 


48728.5 


117.7519 


± 0.0040 


CGRO/BATSE 22 


48733.5 


117.7791 


± 0.0041 


CGRO/BATSE 22 


48738.5 


117.8007 


± 0.0045 


CGRO/BATSE 22 


48743.5 


117.8245 


± 0.0059 


CGRO/BATSE 22 


48748.5 


117.8547 


± 0.0026 


CGRO/BATSE 22 


48758.5 


117.9050 


± 0.0038 


CGRO/BATSE 22 


48763.5 


117.9317 


± 0.0029 


CGRO/BATSE 22 


48768.5 


117.9654 


± 0.0020 


CGRO/BATSE 22 


48773.5 


117.9887 


± 0.0024 


CGRO/BATSE 22 


48778.5 


118.0121 


± 0.0043 


CGRO/BATSE 22 


48788.5 


118.0692 


± 0.0045 


CGRO/BATSE 22 


48793.5 


118.0842 


± 0.0048 


CGRO/BATSE 22 


48808.5 


118.1588 


± 0.0053 


CGRO/BATSE 22 


48813.5 


118.1920 


± 0.0067 


CGRO/BATSE 22 


48818.5 


118.2127 


± 0.0047 


CGRO/BATSE 22 


48823.5 


118.2317 


± 0.0042 


CGRO/BATSE 22 


48828.5 


118.2550 


± 0.0061 


CGRO/BATSE 22 


48833.5 


118.2585 


± 0.0083 


CGRO/BATSE 22 


48838.5 


118.2925 


± 0.0050 


CGRO/BATSE 22 


48843.5 


118.3115 


± 0.0042 


CGRO/BATSE 22 
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Table B.l. Pulse period measurements of GX 1+4 (continued). 



MJD 




P spin (s) 








Instrument 


MJD 




P 'spin (s) 








Instrument 


48848 


.5 


118 


.3342 


± 





.0056 


CGRO/BATSE 22 


48858 


.5 


118 


.3761 


+ 





,0037 


CGRO/BATSE 22 


48863 


.5 


118 


4001 


+ 





.0072 


CGRO/BATSE 22 


48868 


.5 


118 


.4083 


+ 





,0043 


CGRO/BATSE 22 


48873 


,5 


118 


4485 


+ 





.0025 


CGRO/BATSE 22 


48878 


.5 


118 


.4775 


+ 





,0050 


CGRO/BATSE 22 


48883 


,5 


118 


.5092 


+ 


0. 


.0072 


CGRO/BATSE 22 


48888 


.5 


118 


.5331 


+ 





,0045 


CGRO/BATSE 22 


48893 


.5 


118 


.5587 


+ 





.0033 


CGRO/BATSE 22 


48898 


.5 


118 


.5938 


+ 





,0024 


CGRO/BATSE 22 


48903 


,5 


118 


6269 


+ 





.0021 


CGRO/BATSE 22 


48908 


.5 


118 


.6541 


+ 





,0024 


CGRO/BATSE 22 


48913 


,5 


118 


.6938 


± 


0. 


.0024 


CGRO/BATSE 22 


48918 


.5 


118 


.7237 


+ 





,0023 


CGRO/BATSE 22 


48923 


.5 


118 


.7449 


+ 





.0040 


CGRO/BATSE 22 


48928 


.5 


118 


,7706 


+ 





,0037 


CGRO/BATSE 22 


48933 


,5 


118 


.7982 


+ 





.0022 


CGRO/BATSE 22 


48938 


.5 


118 


.8347 


+ 





,0020 


CGRO/BATSE 22 


48943 


,5 


118 


.8681 


+ 


0. 


.0020 


CGRO/BATSE 22 


48948 


.5 


118 


.9047 


+ 





,0020 


CGRO/BATSE 22 


48953 


.5 


118 


.9359 


+ 





.0027 


CGRO/BATSE 22 


48958 


.5 


118 


9641 


+ 





,0045 


CGRO/BATSE 22 


48963 


,5 


119 


0017 


+ 





.0026 


CGRO/BATSE 22 


48968 


.5 


119 


0269 


+ 





,0038 


CGRO/BATSE 22 


48973 


,5 


119 


.0432 


+ 


0. 


0039 


CGRO/BATSE 22 


48978 


.5 


119 


.0783 


+ 





,0029 


CGRO/BATSE 22 


48983 


.5 


119 


1106 


+ 





.0023 


CGRO/BATSE 22 


48988 


.5 


119 


1483 


+ 





0019 


CGRO/BATSE 22 


48993 


,5 


119 


1940 


+ 





.0020 


CGRO/BATSE 22 


48998 


.5 


119 


.2373 


+ 





0019 


CGRO/BATSE 22 


49003 


,5 


119 


.2769 


± 


0. 


.0028 


CGRO/BATSE 22 


49008 


.5 


119 


.3053 


+ 





,0041 


CGRO/BATSE 22 


49013 


.5 


119 


.3430 


+ 





.0039 


CGRO/BATSE 22 


49018 


.5 


119 


.3886 


+ 





,0025 


CGRO/BATSE 22 


49023 


,5 


119 


.4221 


+ 





.0017 


CGRO/BATSE 22 


49028 


.5 


119 


.4549 


+ 





0019 


CGRO/BATSE 22 


49033 


,5 


119 


.4886 


+ 


0. 


.0023 


CGRO/BATSE 22 


49038 


.5 


119 


.5085 


+ 





,0028 


CGRO/BATSE 22 


49043 


.5 


119 


.5367 


+ 





.0044 


CGRO/BATSE 22 


49048 


.5 


119 


.5562 


+ 





,0027 


CGRO/BATSE 22 


49053 


,5 


119 


.5777 


+ 





.0044 


CGRO/BATSE 22 


49058 


.5 


119 


.5973 


+ 





,0041 


CGRO/BATSE 22 


49063 


.5 


119 


6168 


+ 


0. 


.0029 


CGRO/BATSE 22 


49068 


.5 


119 


.6441 


+ 





,0027 


CGRO/BATSE 22 


49073 


.5 


119 


6639 


+ 





.0030 


CGRO/BATSE 22 


49078 


,5 


119 


.6938 


+ 





,0031 


CGRO/BATSE 22 


49083 


,5 


119 


.7111 


+ 





.0047 


CGRO/BATSE 22 


49088 


.5 


119 


.7352 


+ 





,0044 


CGRO/BATSE 22 


49093 


,5 


119 


.7539 


+ 


0. 


.0046 


CGRO/BATSE 22 


49098 


.5 


119 


.7773 


+ 





,0041 


CGRO/BATSE 22 


49103 


.5 


119 


.7868 


+ 





.0078 


CGRO/BATSE 22 


49108 


,5 


119 


,8148 


+ 





,0029 


CGRO/BATSE 22 


49118 


,5 


119 


.8550 


+ 


0. 


.0048 


CGRO/BATSE 22 


49123 


.5 


119 


.8667 


+ 





,0059 


CGRO/BATSE 22 


49128 


,5 


119 


.8897 


+ 


0. 


.0074 


CGRO/BATSE 22 


49133 


.5 


119 


9099 


+ 





0069 


CGRO/BATSE 22 


49138 


.5 


119 


.9292 


+ 





.0039 


CGRO/BATSE 22 


49148 


,5 


119 


,9736 


+ 





,0038 


CGRO/BATSE 22 


49153 


,5 


119 


9993 


+ 





.0055 


CGRO/BATSE 22 


49158 


.5 


120 


0206 


+ 





,0035 


CGRO/BATSE 22 


49163 


,5 


120 


.0446 


+ 


0. 


.0026 


CGRO/BATSE 22 


49168 


.5 


120 


,0767 


+ 





,0018 


CGRO/BATSE 22 


49173 


.5 


120 


.1002 


+ 





.0014 


CGRO/BATSE 22 


49178 


,5 


120 


,1274 


+ 





,0015 


CGRO/BATSE 22 


49183 


,5 


120 


.1575 


+ 





0017 


CGRO/BATSE 22 


49188 


.5 


120 


1903 


+ 





,0018 


CGRO/BATSE 22 


49193 


.5 


120 


.2209 


+ 





0021 


CGRO/BATSE 22 


49198 


.5 


120 


,2545 


+ 





,0027 


CGRO/BATSE 22 


49203 


.5 


120 


3001 


+ 





.0042 


CGRO/BATSE 22 


49208 


.5 


120 


,3250 


+ 





,0057 


CGRO/BATSE 22 


49218 


,5 


120 


.4089 


+ 





.0028 


CGRO/BATSE 22 


49223 


.5 


120 


,4351 


+ 





,0026 


CGRO/BATSE 22 


49228 


.5 


120 


.4876 


+ 





0016 


CGRO/BATSE 22 


49232 


.0 


120 


,5560 


+ 





,0240 


GRANAT/SIGMA 


49233 


.0 


120 


.5570 


+ 





.0280 


GRANAT/SIGMA 23 


49233 


,5 


120 


,5371 


+ 





,0014 


CGRO/BATSE 22 


49236 


,0 


120 


.5730 


+ 





.0150 


GRANAT/SIGMA 23 


49238 


.0 


120 


,5870 


+ 





0150 


GRANAT/SIGMA 


49238 


.5 


120 


.5885 


+ 





0016 


CGRO/BATSE 22 


49242 


.0 


120 


,6440 


+ 





0150 


GRANAT/SIGMA 


49243 


.5 


120 


.6464 


+ 





.0014 


CGRO/BATSE 22 


49244 


.0 


120 


6680 


+ 





,0200 


GRANAT/SIGMA 


49246 


,0 


120 


.6750 


+ 





0080 


GRANAT/SIGMA 23 


49248 


.5 


120 


,7016 


+ 





0030 


CGRO/BATSE 22 


49249 


.0 


120 


.7080 


+ 





0160 


GRANAT/SIGMA 23 


49253 


.5 


120 


,7252 


+ 





,0043 


CGRO/BATSE 22 


49258 


.5 


120 


.7591 


+ 





.0047 


CGRO/BATSE 22 


49263 


,5 


120 


,7738 


+ 





,0035 


CGRO/BATSE 22 


49268 


,5 


120 


.8137 


+ 





.0052 


CGRO/BATSE 22 


49273 


.5 


120 


,8374 


+ 





,0076 


CGRO/BATSE 22 


49288 


.5 


120 


.8429 


+ 





.0043 


CGRO/BATSE 22 


49293 


.5 


120 


,8426 


+ 





,0033 


CGRO/BATSE 22 


49298 


.5 


120 


.8701 


+ 





.0042 


CGRO/BATSE 22 


49303 


,5 


120 


,9023 


+ 





,0045 


CGRO/BATSE 22 


49308 


,5 


120 


.9368 


+ 





.0032 


CGRO/BATSE 22 


49313 


.5 


120 


9719 


+ 





,0024 


CGRO/BATSE 22 


49318 


.5 


121 


.0084 


+ 





.0020 


CGRO/BATSE 22 


49323 


.5 


121 


,0422 


+ 





0019 


CGRO/BATSE 22 


49328 


.5 


121 


.0770 


+ 





.0021 


CGRO/BATSE 22 


49333 


,5 


121 


.1123 


+ 





,0024 


CGRO/BATSE 22 


49338 


.5 


121 


.1465 


+ 





.0021 


CGRO/BATSE 22 


49343 


.5 


121 


1893 


+ 





,0048 


CGRO/BATSE 22 


49363 


.5 


121 


.3394 


+ 





.0043 


CGRO/BATSE 22 


49368 


.5 


121 


,3456 


+ 





,0072 


CGRO/BATSE 22 


49388 


.5 


121 


.3486 


+ 





.0048 


CGRO/BATSE 22 


49398 


.5 


121 


,3689 


+ 





0061 


CGRO/BATSE 22 


49403 


,5 


121 


.3977 


+ 


0. 


.0083 


CGRO/BATSE 22 


49408 


.5 


121 


4013 


+ 





,0070 


CGRO/BATSE 22 


49413 


.5 


121 


.4083 


+ 





.0055 


CGRO/BATSE 22 


49418 


.5 


121 


,4174 


+ 





,0045 


CGRO/BATSE 22 


49423 


.5 


121 


.4414 


+ 





.0030 


CGRO/BATSE 22 


49438 


.5 


121 


,4471 


+ 





,0079 


CGRO/BATSE 22 


49443 


,5 


121 


.4712 


+ 





.0033 


CGRO/BATSE 22 


49448 


.5 


121 


,4846 


+ 





,0026 


CGRO/BATSE 22 


49453 


.5 


121 


.4922 


+ 





.0027 


CGRO/BATSE 22 


49458 


.5 


121 


,5273 


+ 





,0055 


CGRO/BATSE 22 


49468 


.5 


121 


.5630 


+ 





.0073 


CGRO/BATSE 22 


49473 


.5 


121 


,5816 


+ 





,0038 


CGRO/BATSE 22 


49478 


,5 


121 


6012 


+ 





.0033 


CGRO/BATSE 22 


49483 


.5 


121 


,6208 


+ 





0039 


CGRO/BATSE 22 


49488 


.5 


121 


.6330 


+ 





.0045 


CGRO/BATSE 22 


49493 


.5 


121 


,6247 


+ 





,0074 


CGRO/BATSE 22 


49498 


,5 


121 


6371 


± 


0. 


0041 


CGRO/BATSE 22 


49503 


.5 


121 


6693 


+ 





,0062 


CGRO/BATSE 22 


49508 


,5 


121 


.6771 


+ 





.0042 


CGRO/BATSE 22 


49513 


.5 


121 


.7119 


+ 





,0074 


CGRO/BATSE 22 


49518 


.5 


121 


.7488 


+ 





.0057 


CGRO/BATSE 22 


49523 


.5 


121 


,7540 


+ 





,0043 


CGRO/BATSE 22 


49528 


,5 


121 


.7890 


+ 





.0038 


CGRO/BATSE 22 


49533 


.5 


121 


,8093 


+ 





,0038 


CGRO/BATSE 22 
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Table B.l. Pulse period measurements of GX 1+4 (continued). 



MJD 




P spin (s) 








Instrument 


MJD 




P spin (s) 








Instrument 


49538 


.5 


121 


.8352 


+ 





.0029 


CGRO/BATSE 22 


49548 


.5 


121 


.8701 


+ 





0030 


CGRO/BATSE 22 


49558 


.5 


121 


.9098 


+ 





.0022 


CGRO/BATSE 22 


49563 


.5 


121 


.9172 


+ 





0016 


CGRO/BATSE 22 


49573 


,5 


121 


.9533 


+ 





.0027 


CGRO/BATSE 22 


49578 


.5 


121 


9681 


+ 





0036 


CGRO/BATSE 22 


49583 


,5 


121 


.9783 


+ 


0. 


.0029 


CGRO/BATSE 22 


49588 


.5 


121 


9839 


+ 





,0031 


CGRO/BATSE 22 


49593 


.5 


121 


9993 


+ 





0046 


CGRO/BATSE 22 


49598 


.5 


122 


.0037 


+ 





,0033 


CGRO/BATSE 22 


49603 


,5 


122 


0168 


+ 





.0034 


CGRO/BATSE 22 


49608 


.5 


122 


.0263 


+ 





,0023 


CGRO/BATSE 22 


49613 


,5 


122 


.0425 


+ 


0. 


.0025 


CGRO/BATSE 22 


49618 


.5 


122 


0517 


+ 





0019 


CGRO/BATSE 22 


49623 


.5 


122 


.0657 


+ 





.0018 


CGRO/BATSE 22 


49628 


.5 


122 


,0727 


+ 





,0025 


CGRO/BATSE 22 


49633 


,5 


122 


0869 


+ 





0019 


CGRO/BATSE 22 


49638 


.5 


122 


0986 


+ 





,0014 


CGRO/BATSE 22 


49643 


,5 


122 


1063 


± 


0. 


.0014 


CGRO/BATSE 22 


49648 


.5 


122 


1128 


+ 





,0014 


CGRO/BATSE 22 


49653 


.5 


122 


1176 


+ 





.0014 


CGRO/BATSE 22 


49658 


.5 


122 


1159 


+ 





,0012 


CGRO/BATSE 22 


49663 


,5 


122 


1039 


+ 





.0012 


CGRO/BATSE 22 


49668 


.5 


122 


.0851 


+ 





,0012 


CGRO/BATSE 22 


49673 


,5 


122 


0689 


+ 


0. 


.0012 


CGRO/BATSE 22 


49678 


.5 


122 


.0477 


+ 





,0011 


CGRO/BATSE 22 


49683 


.5 


122 


0290 


+ 





0010 


CGRO/BATSE 22 


49688 


.5 


122 


.0089 


+ 





0010 


CGRO/BATSE 22 


49693 


,5 


121 


.9882 


+ 





0011 


CGRO/BATSE 22 


49698 


.5 


121 


.9674 


+ 





,0013 


CGRO/BATSE 22 


49703 


,5 


121 


9506 


+ 


0. 


.0012 


CGRO/BATSE 22 


49708 


.5 


121 


.9345 


+ 





,0012 


CGRO/BATSE 22 


49713 


.5 


121 


.9157 


+ 





.0012 


CGRO/BATSE 22 


49718 


.5 


121 


.8965 


+ 





,0012 


CGRO/BATSE 22 


49723 


,5 


121 


.8798 


+ 





.0012 


CGRO/BATSE 22 


49728 


.5 


121 


.8653 


+ 





,0011 


CGRO/BATSE 22 


49733 


,5 


121 


.8524 


+ 


0. 


0010 


CGRO/BATSE 22 


49738 


.5 


121 


.8382 


+ 





0009 


CGRO/BATSE 22 


49743 


.5 


121 


.8264 


+ 





.0013 


CGRO/BATSE 22 


49748 


.5 


121 


8116 


+ 





0013 


CGRO/BATSE 22 


49753 


,5 


121 


.8013 


+ 





.0014 


CGRO/BATSE 22 


49758 


.5 


121 


.7929 


+ 





,0015 


CGRO/BATSE 22 


49763 


,5 


121 


.7853 


± 


0. 


.0015 


CGRO/BATSE 22 


49768 


.5 


121 


.7783 


+ 





,0011 


CGRO/BATSE 22 


49773 


.5 


121 


.7763 


+ 





.0011 


CGRO/BATSE 22 


49778 


.5 


121 


.7751 


+ 





,0015 


CGRO/BATSE 22 


49783 


,5 


121 


.7749 


+ 





.0012 


CGRO/BATSE 22 


49788 


.5 


121 


.7792 


+ 





,0012 


CGRO/BATSE 22 


49793 


,5 


121 


.7886 


+ 


0. 


.0013 


CGRO/BATSE 22 


49798 


.5 


121 


.7998 


+ 





,0017 


CGRO/BATSE 22 


49803 


.5 


121 


.8107 


+ 





0019 


CGRO/BATSE 22 


49808 


.5 


121 


.8269 


+ 





,0028 


CGRO/BATSE 22 


49813 


,5 


121 


.8397 


+ 





.0045 


CGRO/BATSE 22 


49818 


.5 


121 


.8633 


+ 





0069 


CGRO/BATSE 22 


49823 


,5 


121 


.9408 


± 


0. 


0066 


CGRO/BATSE 22 


49833 


.5 


122 


.0374 


+ 





,0024 


CGRO/BATSE 22 


49838 


.5 


122 


1106 


+ 





.0017 


CGRO/BATSE 22 


49843 


,5 


122 


1481 


+ 





,0045 


CGRO/BATSE 22 


49848 


,5 


122 


1751 


+ 





.0020 


CGRO/BATSE 22 


49853 


.5 


122 


.2054 


+ 





0019 


CGRO/BATSE 22 


49858 


,5 


122 


.2373 


± 


0. 


.0017 


CGRO/BATSE 22 


49863 


.5 


122 


.2816 


+ 





,0021 


CGRO/BATSE 22 


49868 


.5 


122 


.3200 


+ 





.0017 


CGRO/BATSE 22 


49873 


,5 


122 


,3576 


+ 





,0030 


CGRO/BATSE 22 


49878 


,5 


122 


3916 


+ 





.0029 


CGRO/BATSE 22 


49883 


.5 


122. 


.4229 


+ 





,0017 


CGRO/BATSE 22 


49888 


.5 


122 


.4551 


+ 





.0015 


CGRO/BATSE 22 


49893 


.5 


122 


,4858 


+ 





,0017 


CGRO/BATSE 22 


49898 


.5 


122 


.5197 


+ 





.0021 


CGRO/BATSE 22 


49903 


.5 


122 


,5389 


+ 





,0022 


CGRO/BATSE 22 


49908 


,5 


122 


.5647 


+ 





0016 


CGRO/BATSE 22 


49913 


.5 


122. 


,5870 


+ 





,0017 


CGRO/BATSE 22 


49918 


.5 


122 


6001 


+ 





.0018 


CGRO/BATSE 22 


49923 


.5 


122 


,6279 


+ 





0016 


CGRO/BATSE 22 


49928 


.5 


122 


.6701 


+ 





0019 


CGRO/BATSE 22 


49933 


,5 


122 


6951 


+ 





,0014 


CGRO/BATSE 22 


49938 


,5 


122 


.7172 


+ 





0016 


CGRO/BATSE 22 


49943 


.5 


122. 


,7428 


+ 





,0022 


CGRO/BATSE 22 


49948 


.5 


122 


.7698 


+ 





0036 


CGRO/BATSE 22 


49953 


.5 


122 


,7684 


+ 





,0015 


CGRO/BATSE 22 


49958 


,5 


122 


.7916 


+ 





0016 


CGRO/BATSE 22 


49963 


,5 


122 


,8126 


+ 





,0014 


CGRO/BATSE 22 


49968 


,5 


122 


.8554 


+ 





0019 


CGRO/BATSE 22 


49971 


.0 


122 


,8770 


+ 





,0210 


GRANAT/SIGMA 


49973 


.5 


122 


.8844 


+ 





0021 


CGRO/BATSE 22 


49978 


.0 


122 


,9230 


+ 





0090 


GRANAT/SIGMA 


49978 


.5 


122 


.9032 


+ 





.0013 


CGRO/BATSE 22 


49979 


.0 


122 


,9420 


+ 





,0400 


GRANAT/SIGMA 


49983 


,5 


122 


.9475 


+ 





.0015 


CGRO/BATSE 22 


49988 


.5 


122 


,9744 


+ 





,0017 


CGRO/BATSE 22 


49993 


.5 


122 


9989 


+ 





0016 


CGRO/BATSE 22 


49998 


.5 


123 


0159 


+ 





,0015 


CGRO/BATSE 22 


50003 


.5 


123 


.0368 


+ 





.0013 


CGRO/BATSE 22 
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127. 


.7850 


+ 





.0015 


CGRO/BATSE 29 


50959 


.9 


127. 


,7770 ± 


0.0015 


CGRO/BATSE 29 


50964 


1 


127. 


8110 


± 





.0015 


CGRO/BATSE 29 


50968 


.0 


127. 


,8420 ± 


0.0015 


CGRO/BATSE 29 


50972 


.0 


127. 


.8710 


± 





0011 


CGRO/BATSE 29 


50976 


.0 


127. 


,9080 ± 


0.0011 


CGRO/BATSE 29 


50980 


,0 


127. 


.9390 


± 


0. 


.0007 


CGRO/BATSE 29 


50984 


.0 


127. 


,9730 ± 


0.0007 


CGRO/BATSE 29 


50988 


,0 


128 


.0080 


± 





0006 


CGRO/BATSE 29 


50992 


.1 


128 


,0420 ± 


0.0006 


CGRO/BATSE 29 


50996 


.0 


128 


.0750 


+ 





.0007 


CGRO/BATSE 29 


50999 


.5 


128 


1000 ± 


0.0007 


CGRO/BATSE 29 


51004 


.8 


128 


1520 


+ 





.0015 


CGRO/BATSE 29 


51008 


,0 


128 


1780 ± 


0.0015 


CGRO/BATSE 29 
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Table B.l. Pulse period measurements of GX 1+4 (continued). 



MJD 




P spin (s) 






Tnstriimpnt 




MLD 




spinas) 






Tir;ti"iimpnt 

Ll Lll 1 IV. 1 J I. 


51012 




128 


.2140 ± 





.0007 


CGRO/RATSF ' a 

V V J J.VVJ/ LJI\ 1 k ) Ij 




51016. 


,0 


128 


,2520 ± 





,0007 


CGRO/RATSF 29 

VAI1VV// IJ J \ 1 Jtj 


51020 


1 


128 


.2880 ± 


n 


0007 
uuu i 


PPRPj/RATSF 29 

V^VJIvVJ/ IJ rt 1 rtLrt 




51073 

J 1UZJ 


q 


128 


,3210 ± 


A 
yt 


0007 
uuu / 


PPtRO/RATSF 29 

V . V ) 1 VV )l LJ rt 1 OLrt 


51028 


.0 


128 


.3560 + 





.0007 


CGRO/RATSF 29 

V V J IX V 7/ LJ I \ 1 k) Ij 




51032. 


o 


128 


,3920 ± 





,0007 


CGRO/RATSF 29 

V. VJ1VV7/ IJ ^ \ 1 JLi 


51036 


.0 


128 


.4270 ± 


o 


0006 


CGRO/RATSF 29 

V V J IX V Si LJ I \ 1 k) I j 




51040. 


o 


128 


,4640 ± 





0006 


CGRO/RATSF 29 

v-vjivw/urii oij 


51044 





128 


.5010 ± 


n 


0004 
yjyjyj L + 


PPRO/RATSF 29 

V rt J IV v 7/ LJ rt 1 rtLrt 




51046 




128 


,5250 ± 


A 
yt 


0004 

UUUt - 


PPtRO/RATSF 29 

VAI 1 vv )l lj rt 1 OLrt 


J 1 U^T / 


c 
o 


128 


.5442 ± 


Q 


0107 


RpnnnSAY/WFP 25 

DCUUUurtA/ VV I V. 




5105? 


q 

7 


128 


,5840 ± 


o 


0008 


PPtRO/RATSF 29 

V . V I L vv si LJ rt L OJTj 


5105IS 




128 


.6140 + 


n 
\t 


0008 
uuuo 


PPtRPj/RATSF 29 

V A l IVA )l LJ rt 1 O-Ci 




51057 

J 1UJ / 


8 
o 


128 


,6338 ± 


A 
u 


0053 


RpnnnSAYAVFP 25 

D CUUUOrtA/ VV I V 


5 1 060 

J 1UUU 


1 


128 


.6530 ± 


A 

U, 


OOOS 


mpo/RAT^F 29 




S 1 On 1 * 


q 

7 


128 


,6920 ± 


A 

yt 


000 S 

UUU J 


PPRO/RAT^F 29 

V rt ) Ivv si Drt 1 rt Irt 


5 1 067 

J 1UU / 


z 


128 


.7238 ± 


A 
U 


004q 

UU47 


Rpnnn>AYAA/Ff 25 




S 1 0n8 
J 1UOO 


A 
U 


128 


,7300 ± 


A 

u 


0004 

UUUH- 


PPRO/RATSF 29 

V . V J IVV SI Drt 1 rtLrt 


5107? 

J 1U / Z 


J 


128 


.7680 ± 


A 
u 


0004 


rriRO/RATSF 29 




51075 


q 


128 


,8050 ± 


A 
yt 


0004 

UUUt - 


PPtRO/RATSF 29 

V rt I IVV SI LJ rt 1 rtJrt 


5 1 07Q 


s 
o 


128 


.8410 ± 


A 

U, 


0004 


phro/rat^f 29 




S 1 07q 




z 


128 


,8189 ± 


A 
u 


01 77 

U1Z / 


RpnnnSAYAVFf 25 

DCppUrt rtyv/ Wrv. 


S 1 084. 


a 

V/ 


128 


.8820 ± 


A 

U, 


0003 


PHRD/RAT^F 29 




S 1 084 

J 1UOH 


u 


128 


,8474 ± 


A 
U 


0740 

UZ4U 


Rpnnn>AYAA/Fr 25 


S 1 088 

J 1UOO 


a 

U 


128 


.9190 ± 


A 

U, 


0003 
uuuj 


phro/rat^f 29 




S 1 00? 

j i yiy z.. 


A 
V7 


128 


,9570 ± 


A 

yt 


000^ 
uuu J 


PPRO/RAT^F 29 

v rt J Ivv si Drt 1 rtJrt 


S 1 OQ? 


J 


128 


.9581 + 


A 

u, 


0033 


D pnn ,,e a V /WFC" 25 
Dcppuort-A./ WrL 




^ 1 0q6 

2)1 uyo 


A 

u 


128 


,9920 ± 


A 

u, 


OOO^ 


rr.Rfl/RAT^F 29 
vrt i IvA 7*/ Drt 1 rtJrt 


S 1 OQQ 

J 1U77 


Q 


129 


.0240 + 


A 
U 


oon? 
uuuz 


rr.pn/RATQP 29 

vAi r\V j/ on 




^110^ 


Q 
o 


129 


,0580 ± 


A 

yt 


000? 
uuuz 


rr,PD/RAT<vF 29 
vrt i IvA L/ Drt 1 rtJrt 


S 1 1 07 

j 1 1 u / 


a 


129 


.0890 + 


A 

U, 


OOOn 

uuuu 


PHRO/RAT^F 29 




S 1 1 1 ? 


A 
V7 


129 


.1230 + 


A 

yt 


OOOn 
UUUO 


PPRO/RAT^F 29 

V rt J Ivv si Drt 1 Dirt 


S 1 1 1 f\ 


A 
17 


129 


.1520 + 


A 
U 


0008 
UUUO 


rr,pn/RAT<vP 29 




S 1 1 ?0 

2> 1 IZu. 


1 
1 


129 


,1840 ± 


A 

u, 


0008 
UUUo 


rr.Rfl/RAT^F 29 
vrt 1 IvV L/ Drt 1 rtLrt 


S 1 1 74 


A 
17 


129 


.2150 ± 


A 
U 


001 
UU1U 


rr T i?n/RAT<vP 29 

vAi r\V J/ oil 




^1178 
2> 1 IZo 


1 
1 


129 


,2460 ± 


A 

yt 


001 
UU 1U 


rr,pn/RAT<vF 29 

vrt ) IvA L/ Drt 1 oJrt 


S 1 1 3? 


A 
At 


129 


.2780 ± 


A 

U, 


0008 

UUUO 


PPtRO/RAT^F 29 

vAI I\A ^/ Drt 1 rt Irt 




SI 1 3fi 


1 
1 


129 


3100 ± 


A 

yt 


0008 

UUUO 


PPRO/RAT^F 29 

V rt J Ivv si Drt 1 rt Irt 


51 1 39 


q 

y 


129 


.3420 ± 


A 
u 


0008 

UUUO 


rHRO/RATSF 29 

V A I I XV )l LJ r\ 1 O-Ci 




51 144 




129 


,3740 ± 


A 
yt 


0008 

UUUO 


PPtRO/RATSF 29 

V . V J IVV SI LJ rt 1 rt Jrt 


S 1 1 47 

J 1 14/ 


q 

7 


129 


4040 ± 


A 

U, 


0008 

UUUO 


PPtRO/RAT^F 29 




S1 1 5? 

J 1 1 JZ 


A 
At 


129 


,4420 ± 


A 
yt 


0008 

UUUO 


PPRO/RAT^F 29 

vrt) Ivv si Drt 1 rt Irt 


S 1 1 Sn 

J 1 1 JU 


A 
At 


129 


.4770 ± 


A 

U, 


OOOS 

UUUJ 


PPtRO/RAT^F 29 

v rt J I\v Drt 1 rt Irt 




S 1 1 nO 

J 1 1 uu 


A 
At 


129 


.5110 ± 


A 

yt 


000 S 

UUU J 


PPRO/RAT^F 29 

V rt J Ivv si Drt 1 rt Irt 


51 1n3 


e 
o 


129 


.5450 ± 


A 
u 


OOOQ. 

UUU7 


PPRO/RATSF 29 

vAI I XV )l LJ rt 1 O-Ci 




51 1(S8 

■J 1 1 UO 


A 
U 


129 


,5790 ± 


A 
yt 


000Q 
yjyjyjy 


PPtRO/RATSF 29 

V . V J IVV SI LJ rt 1 rtLrt 


S 1 1 7? 


A 
At 


129 


6100 ± 


A 

U, 


0008 

UUUO 


PPRn/RAT*\F 29 

v rt J L\v 7/ Drt 1 ljJ3 




S 1 1 7n 

J 1 1 / u 


A 
At 


129 


,6470 ± 


A 
yt 


0008 

UUUO 


PPRO/RAT^F 29 

vrt) Ivv si Drt 1 rt Irt 


51180 


.0 


129 


.6830 ± 


o 


.0008 


PPtRO/RATSF 29 

V A 1 1\V7/ LJ / \ 1 k ) I 




51 184. 





129 


,7190 ± 





,0008 


PGRO/RATSF 29 

V . V J IVV SI IJ J \ 1 JLi 


51 1 88 


1 


129 


.7520 ± 


A 


0007 
uuu / 


PPRO/RATSF 29 

V VI I VV )l LJ rt 1 O-Cj 




51 1Q? 

Jll 7i. 


A 
U 


129 


,7870 ± 


A 
u 


0007 
uuu / 


PPtRO/RATSF 29 

V rt J 1 VV SI LJ rt 1 rtLrt 


51196 


1 


129 


.8230 ± 


o 


.0008 


TGRO/RATSF 29 




51200 





129 


,8580 ± 





,0008 


CGRO/RATSF 29 

V. VJ1VV7/ IJ f \ 1 JLj 


51204 


o 


129 


.8940 ± 





.0005 


CGRO/RATSF 29 




51208. 


o 


129 


9310 ± 





,0005 


CGRO/RATSF 29 


51212 


.0 


129 


.9650 ± 





.0005 


CGRO/RATSF 29 




51216. 


o 


130 


,0030 ± 





,0005 


CGRO/RATSF 29 

V.A J IW7/ LJ J \ 1 Jlj 


5171Q 


q 


130 


.0370 ± 


A 
u 


0005 


mRO/RATSF 29 




51224 


A 

At 


130 


,0750 ± 


o 


0005 


CGRO/RATSF 29 

V 1, )I\V !) \ 1 ) 1 . 


51228 


q 


130 


,1120 ± 


A 

U, 


0005 

UUUJ 


rnRn/RATSF 29 




51730 

■J 1ZJU 


3 


130 


,1523 ± 


A 
yt 


0074 


RpnnnSAY/WFC 25 

iJCULJUOrtTV/ vv rv. 


5173? 

J 1ZJZ 


q 


130 


.1530 + 


Q 


0005 


mRO/RATSF 29 




51?3n 


1 


130 


,1930 ± 


o 


0005 


CGRO/RATSF 29 


51240 


q 


130 


.2320 ± 


Q 


0005 


mRD/RATSF 29 




51244 


1 


130 


,2730 ± 


o 


0004 


CGRO/RATSF 29 


51748 


1 


130 


.3150 + 


A 

U, 


0004 


rnRn/RATSF 29 




51740 


4 


130 


,3295 ± 


A 
yt 


0140 

UltU 


RpnnnSAY/WFC 25 

iJCULJUOrtTV/ vv rv. 


51751 


q 


130 


.3560 ± 


A 
u 


0004 


mRO/RATSF 29 




51?5n 


A 
U 


130 


,3990 ± 


Q 


0004 


CGRO/RATSF 29 


S 1 7SQ 


q 

7 


130 


.4410 ± 


A 
u 


0003 

UUUJ 


PHRD/RATSF 29 




S 1 ?n? 
J 1ZUZ 


J 


130 


,4689 ± 


A 
u 


01 ^1 


RpnnnSAX/WFC 25 


S 1 7n3 


q 

7 


130 


4880 ± 


A 

U, 


0003 

UUUJ 


PHRD/RAT^F 29 




S 1 7n7 

J 1 zu / 


q 

7 


130 


,5330 ± 


A 

yt 


000^ 

UUUJ 


CGRO/RATSF 29 


51270 


8 


130 


.5625 ± 





.0074 


RennoSAXAVFC 25 




51272. 





130 


,5780 ± 





,0003 


CGRO/RATSF 29 

V.A J l\V// LJ I V 1 i 1 1 j 


51274 


7 


130 


.6058 ± 


Q 


OOQQ 
yjytyy 


RpnnnSAX/WFP 25 
d eppuo rt.yv/ vv iv. 




51?7n 




130 


,6250 ± 





0003 


PPRO/RATSF 29 

V rt 1 Ivv )l LJ rt 1 rtLrt 


S 1 778 


.j 


130 


.6510 + 


A 
u 


0047 

AJU4 / 


Rpnnn^AYAVFf 25 
ocppuortyv/ vvrv. 




S 1 780 

J 1ZOU 


1 
1 


130 


,6720 ± 


A 

u 


0003 

UUUJ 


PPRO/RATSF 29 

vA) Ivv )l Drt 1 rt Jrt 


51784 


1 


130 


.7200 ± 


A 
yt 


0003 

UUUJ 


PPRO/RATSF 29 

V rt 1 1 vv )l LJ rt 1 OJ_/ 




51788 

^> 1ZOO 


A 
At 


130 


,7660 ± 


A 

yt 


0003 

UUUJ 


PPtRO/RATSF 29 

V rt 1 1 vv )l IJ rt 1 rtLrt 


51292 




130 


.8130 + 


A 
yt 


000? 
uuuz 


PPRO/RATSF 29 

V rt J 1\V Si LJ rt 1 JEi 




J 1Z7U 




130 


,8590 ± 


() 


0002 


PPRO/RATSF 29 

V rt J Ivv )l LJ rt 1 rtLrt 


51 300 

J UUU 


1 


130 


.9050 ± 


A 
yt 


000? 
uuuz 


PPRO/RATSF 29 

V rt 1 1 vv )l LJ rt 1 O-Ci 




51 304 


A 
At 


130 


,9530 ± 


A 
yt 


0002 


PPtRO/RATSF 29 

V . V j ivv )l Ij rt 1 rtLrt 


51308 


.0 


130 


9990 ± 





.0002 


TGRD/RATSF 29 




51312. 





131 


,0450 ± 





,0002 


CGRO/RATSF 29 

V.A J IVV7/ LJ J \ 1 JLj 


51315 


q 


131 


0910 ± 





.0003 


CGRO/RATSF 29 

V V J 1 V V s/LJ / \ 1 k ) I 




51320. 


o 


131 


,1390 ± 





,0003 


CGRO/RATSF 29 

V. VI1VV7/ IJ ^ \ 1 JLj 


51324 


.0 


131 


.1850 ± 





.0003 


CGRO/RATSF 29 

V viivvy/ij / \ i 




51328 


o 


131 


,2280 ± 





,0003 


CGRO/RATSF 29 

V.A J IVV7/ LJ J \ 1 JLj 


51332 


o 


131 


.2740 ± 





0003 


CGRO/RATSF 29 

vvjivvy/ lj i v i tii i 




51336. 


o 


131 


,3210 ± 





,0003 


CGRO/RATSF 29 

V.A J l\V// LJ I V 1 i J 1 j 


51340 


o 


131 


.3650 ± 





.0004 


CGRO/RATSF 29 

V V J 1 VV St LJ I V 1 l l 1 i 




51344. 


o 


131 


,4090 ± 





,0004 


CGRO/RATSF 29 


51460 


q 


132 


.5723 ± 





0089 


RennoSAXAVFC 25 




51464. 


3 


132 


,6220 ± 





,0174 


RennoSAXAVFC 25 

jU L^ IJ UUvJil/V/ VV LY V 


51785 


o 


134 


.9256 ± 


0. 


0010 


BeppoSAX 24 




52346. 


1 


137 


,4690 ± 





,3680 


RXTF/PCA&HFXTF 26 

IV / V J. J il L V I V VX. 1 1L//V LV 1 i 


52351 


3 


137 


.2540 ± 





2910 


RXTE/PCA&HEXTE 


26 


52359. 


q 


137 


,7750 ± 





,0920 


RXTF/PCA&HFXTF 26 

XVjtV LV J il L V 1 V VV_ 1 1 1 j X L_y 


52366 


.6 


137 


.9440 ± 





2760 


RXTE/PCA&HEXTE 


26 


52384. 


7 


137 


,8050 ± 





1690 


RXTF/PCA&HFXTF 26 

IV / V J. J il L V I V VX. 1 1L//V LV 1 i 


52392 


3 


1 11 


neon _i_ 





1840 


RXTE/PCA&HEXTE 


26 


52401 


2 


Lj l 


OUJU ± 





,2000 


RXTF/PCA&HFXTF 26 


52407 


,0 


137 


.6350 ± 





1840 


RXTE/PCA&HEXTE 


26 


52418. 





138 


0180 ± 





,2610 


RXTE/PCA&HEXTE 26 


52422 


.0 


137 


.1890 + 





2910 


RXTE/PCA&HEXTE 


26 


52431. 


,2 


138 


,2020 ± 


0, 


,4760 


RXTE/PCA&HEXTE 26 


52434 


,9 


138 


0790 ± 





1990 


RXTE/PCA&HEXTE 


26 


52450. 


8 


137, 


,5860 + 


1 


,8740 


RXTE/PCA&HEXTE 26 


52466 


1 


137 


1250 + 





.3070 


RXTE/PCA&HEXTE 


26 


52470. 


,7 


137 


9390 ± 


0, 


,1840 


RXTE/PCA&HEXTE 26 


52490 


,9 


138 


.1370 + 


0. 


.2300 


RXTE/PCA&HEXTE 


26 


52491. 





138 


,1700 + 





0010 


Chandra&RXTE 27 


52493 


.0 


138 


.1680 + 





.1230 


RXTE/PCA&HEXTE 


26 


52499. 


,2 


138 


1060 + 


0, 


1690 


RXTE/PCA&HEXTE 26 


52502 


.2 


138 


.1060 + 


0. 


.5530 


RXTE/PCA&HEXTE 


26 


52503. 


5 


138 


,2590 ± 





,8450 


RXTE/PCA&HEXTE 26 


52512 


.6 


138 


.7040 ± 





4150 


RXTE/PCA&HEXTE 26 


52515, 


,7 


138 


,3350 + 


0, 


,1850 


RXTE/PCA&HEXTE 26 


52521 


,8 


138 


.3660 ± 


0. 


4610 


RXTE/PCA&HEXTE 


26 


52526. 


1 


138 


,3960 ± 





,4920 


RXTE/PCA&HEXTE 26 


52528 


.9 


138 


.0580 ± 





.6760 


RXTE/PCA&HEXTE 


26 


52535, 





138 


1350 + 


0, 


2610 


RXTE/PCA&HEXTE 26 


52539 


,6 


138 


4570 ± 





.3840 


RXTE/PCA&HEXTE 


26 


52542, 


3 


138 


,1650 + 





,6450 


RXTE/PCA&HEXTE 26 


52548 


4 


137 


.9960 ± 





.5680 


RXTE/PCA&HEXTE 


26 


52555, 


,5 


138 


,0570 ± 


3, 


,2090 


RXTE/PCA&HEXTE 26 
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Table B.l. Pulse period measurements of GX 1+4 (continued). 



MJD 




P spin (s) 








Instrument 




MJD 










Instrument 




52561 


o 


138 


.0570 




1 


5660 


R XTF /PC A & H F XTF 


26 


52565 


.6 


138 


.1180 ± 


0. 


7980 


RXTE/PCA&HEXTE lb 


52582 




138 


1630 




o 


.4920 


R XTF/Pr A # HF XTF 

IV /\ J. Ij/ A V /V C>L J 1 1 j/y 1 L j 


26 


52585 


.8 


138 


.6700 ± 


0. 


2910 


RXTE/PCA&HEXTE 26 


52589 


2 


138 


.6230 







1080 


R XTF/Pr A # HF XTF 

IV /\ J. J_,/ L V y l\ OC J 1 Lj/\ 1 Li 


26 


52591 


9 


138 


.5770 ± 


0. 


2000 


RXTE/PCA&HEXTE 26 


52594 


1 


138 


.4850 




o 


.0920 


R XTF/Pr A # HF XTF 


26 


52917 


.6 


139 


.6300 ± 


0. 


0060 


txt 1 i "i 1 n at ^8 

INTEGRAL 




53052 


1 


140 


.6132 


+ 


0. 


.0002 


TNTFGRAT 28 

1.1 v i i j v j i v i v i i 




53252 


4 


141 


.5649 ± 


0. 


0001 


TXT 1 ' 1 .' Vt AT ''S 

IN 1 buKAL 




53422 


o 


141 


3310 


+ 





.3380 


TNTFGRAT /TSGRT 25 

1.1 V 1 l- j V I 1\ I V 1 </ 1 1. J v 1 1\ 1 




53431 


.0 


140 


.9101 ± 


0. 


5285 


TXT 1 ' lj /— " n AT /11VD1 


25 


53434 


o 


141 


.3730 


+ 





.2320 


TNTFGRAT /TSGRT 25 

1.1 i 1 IjVJlvy VI i! 1 k) VI 1\ 1 




53464 


.0 


141 


.8630 ± 


0. 


2986 


TXT 1 ' 1 T1 AT /TC/^FJT 


25 


53467 


.0 


141 


.6670 


+ 


o 


1540 


TNTFGRAT /TSGRT 25 

1.1 1 L J i V I 1 V /V J // IklVIlVl 




53470 


.0 


141 


.4530 ± 


0. 


2489 


TXT 1 ' 1 VI AT /11VD1 

IN IbuKAb/lbuKl 


25 


53473 


o 


141 


.5861 


+ 


o 


.3421 


TNTFGRAT /TSGRT 25 

1.1 1 1 1 jv I IV i V Lj/ 1 k ) V 1 1\ 1 




53476 


.0 


141 


.6270 ± 


0. 


1014 


TXTT'T-7/~ , TJ AT /TCTDT 

IN IbvjKAb/lbCjKl 


25 


53479 


.0 


141 


.8120 


+ 





1796 


TNTFGRAT /TSGRT 25 

1.1 1 L J i V I 1 V /V J // IkJVJIVl 




53611 


.0 


142 


.9630 ± 


0. 


2911 


TXT'T v Lj , f~ , D AT /TCTDT 

IN IbuKAb/lbOKl 


25 


53620 


o 


142 


9510 


+ 


0. 


.3223 


TNTFGRAT /TSGRT 25 

1.1 1 1 IjVIIV/ VLj/ 1 k) VI 1\ 1 




53635 


.0 


142 


.8820 ± 


0. 


1238 


TXTT'T-T/"' D AT /TCTDT 

IN IbuKAb/loOKl 


25 


53638 


.0 


143 


1800 


+ 


o 


1695 


TNTFGRAT /TSGRT 25 

1.1 1 L L j V I IV I V 1 </ 1 1.1 V I IV 1 




53641 


.0 


143 


.0350 ± 


0. 


2270 


IXTTCrD AT /1CVD I 

IN IbuKAb/loOKl 


25 


53644 


,0 


143 


.0580 


+ 


0. 


.2625 


INTEGRAL/ISGRI 25 




53647 


.0 


142 


.6760 ± 


0. 


7286 


TXTTTJ/^O AT /ICY'D I 

IN IbuKAb/lboKl 


25 


53650 


.0 


142 


.9760 


+ 





.3018 


INTEGRAL/IS GRI 25 




53653 


.0 


143 


.0450 ± 


0. 


1820 


TMT'Tir^'O AT /TCTDT 

IN IbuKAb/loOKl 


25 


53656 


,0 


143 


1370 


+ 





.2403 


INTEGRAL/ISGRI 25 




53659 


.0 


143 


1670 ± 


0. 


2502 


TXTTTJ/^LJ AT /TCPDT 

IN IbuKAb/lbOKl 


25 


53662 


.0 


143 


.2370 


+ 





1376 


INTEGRAL/ISGRI 25 




53665 


.0 


143 


.2250 ± 


0. 


2702 


TMTCrD AT /1CVD I 

IN IbuKAb/loOKl 


25 


53668 


,0 


143 


.0780 


+ 


0. 


1641 


INTEGRAL/ISGRI 25 




53776 


.0 


143 


.7540 ± 


0. 


4973 


TXTTTJ/^O AT /ICY'D I 
IN IbvjKAb/lbvjKl 


25 


53779 


.0 


143 


.8440 


+ 





1433 


INTEGRAL/ISGRI 25 




53782 


.0 


143 


,8710 ± 


0. 


1307 


TMTCrD AT /TCTDT 

IN IbuKAb/lbOKl 


25 


53785 


,0 


143 


.8310 


+ 





1705 


INTEGRAL/ISGRI 25 




53788 


.0 


143 


.9010 ± 


0. 


0930 


TXTTTJ/^O AT /TCTDT 

IN IbuKAb/lbLjKl 


25 


53791 


.0 


143 


9800 


+ 





.0914 


INTEGRAL/ISGRI 25 




53794 


.0 


144 


.0310 ± 


0. 


1641 


TMTCrD AT /1CVD I 

IN IbuKAb/loOKl 


25 


53797 


,0 


144 


0950 


+ 


0. 


.2588 


INTEGRAL/ISGRI 25 




53803 


.0 


144 


.0920 ± 


0. 


1138 


TXTTTJ/^O AT /TCTDT 

IN IbuKAb/lboKl 


25 


53806 


.0 


144 


.2470 


+ 





1827 


INTEGRAL/ISGRI 25 




53809 


.0 


144 


.0340 ± 


0. 


1718 


IMTCrD AT /TCTDT 

IN IbuKAb/lbOKl 


25 


53812 


,0 


144 


1330 


+ 





1196 


INTEGRAL/ISGRI 25 




53815 


.0 


144 


.0920 ± 


0. 


1236 


TXTT'T-T/"' D AT /ICY'D I 

IN IbuKAb/lbOKl 


25 


53818 


.0 


144 


.5760 


+ 





.2612 


INTEGRAL/ISGRI 25 




53821 


.0 


144 


.2650 ± 


0. 


1599 


TMTCrD AT /1CVD I 

IN IbuKAb/loOKl 


25 


53830 


.0 


144 


.2780 


+ 





1662 


TNTFGRAT /TSGRT 25 

111 L L j V I IV I V 1 </ 1 i.I V I IV 1 




53833 


.0 


144 


1490 ± 


0. 


1517 


TXT 1 ' 1 VI AT /11VD1 

IN IbOKAb/lbuKl 


25 


53836 


o 


144 


.2310 


+ 


0. 


.2570 


TNTFGRAT /TSGRT 25 

1.1 1 I I j VI 1 V I V t j/ IklVIlVl 




53839 


.0 


144 


.2210 ± 


0. 


1294 


TXT 1 ' 1 .' VI AT / 1 1 • /" T~l I 

IN IbuKAb/lSCjRl 


25 


53842 


.0 


144 


.4540 


+ 





1796 


TNTFGRAT /TSGRT 25 

1.1 1 L I j V 1 1 V I V I _// IklVIlVl 




53845 


.0 


144 


.4040 ± 


0. 


0892 


TXT 1 ' 1 VI AT /11VD1 

IN IbuKAb/lSuKl 


25 


53965 


o 


145 


.3250 


+ 


0. 


.2286 


TNTFGRAT /TSGRT 25 

1.1 1 I I j V 1 1 V I V Lj/ IklVIlVl 




53968 


.0 


145 


.5370 ± 


0. 


1574 


TXT 1 ' 1 .' VI AT / 1 1 • /" ' F) I 

IN IbuKAb/lSCjRl 


25 


53974 


.0 


145 


.4100 


+ 


o 


1504 


TNTFGRAT /TSGRT 25 

1.1 1 L L j V I IV I V 1 </ IklVIlVl 




53977 


.0 


145 


.5250 ± 


0. 


1218 


TXT 1 ' 1 VI AT /11VD1 

IN IbuKAb/lbuKl 


25 


53980 


o 


145 


.4940 


+ 


0. 


1308 


TNTFGRAT /TSGRT 25 

1.1 1 I IjVIiVy V L j/ IklVIlVl 




53983 


.0 


145 


.5390 ± 


0. 


0991 


TXT 1 ' 1 VI AT / 1 1 • /" ' ri I 

IN IbuKAb/lSCjRl 


25 


53986 


.0 


145 


.5590 


+ 





.0952 


TNTFGRAT /TSGRT 25 

1.1 1 L I j V 1 1 V I V 1 </ IklVIlVl 




53989 


.0 


145 


.7140 ± 


0. 


0863 


TXT 1 ' 1 VI AT /11VD1 

IN IbuKAb/lbuKl 


25 


53992 


o 


145 


.7020 


+ 





.0868 


TNTFGRAT /TSGRT 25 

1.1 1 L I j V. 1 1 V I V 1 it I v. I V 1 LV J. 




53998 


.0 


145 


.7880 ± 


0. 


1160 


TATTCr'n AT / 1 1 • /" ' F) I 

IN IbuKAb/lSCjRl 


25 


54001 


.0 


145 


.6570 


+ 





1365 


TNTFGRAT /TSGRT 25 

1.1 1 L L j V I IV I V 1 </ IklVIlVl 




54004 


.0 


145 


,8051 ± 


0. 


1143 


TXT 1 ' 1 ' VI AT /TcrriT 

IN IbuKAb/lSuKl 


25 


54010 


o 


145 


9100 


+ 


0. 


.0883 


TNTFGRAT /TSGRT 25 

1.1 1 I I j VI 1 V i V Lj/ IklVIlVl 




54013 


.0 


145 


.8280 ± 


0. 


1115 


TXT 1 ' 1 VI AT / 1 1 • /" ' ri I 

IN IbuKAb/lSCjRl 


25 


54022 


.0 


145 


.8840 


+ 





.0907 


TNTFGRAT /TSGRT 25 

1.1 1 L 1 j V 1 IV I V I _// IklVIlVl 




54025 


.0 


146 


,0170 ± 


0. 


0638 


TXT 1 ' 1 VI AT /TC/^FJT 

IN IbuKAb/lbuKl 


25 


54028 


o 


146 


1190 


+ 


o 


.0735 


TNTFGRAT /TSGRT 25 

1.1 1 I IjVllVy VI it IklVIlVl 




54031 


.0 


146 


,0260 ± 


0. 


0850 


TXT 1 ' 1 .' VI AT / 1 1 • /" ' F) I 

IN IbuKAb/lSCjRl 


25 


54034 


.0 


146 


0960 


+ 





.0748 


TNTFGRAT /TSGRT 25 

111 L L j V I IV I V 1 </ IklVIlVl 




54148 


.0 


147 


,2450 ± 


0. 


0681 


TXT 1 ' 1 VI AT /I tv D I 

IN IbuKAb/lbuKl 


25 


54151 


o 


147 


.3313 


+ 


0. 


.0736 


TNTFGRAT /TSGRT 25 

1.1 1 I IjVllVy V L j/ IklVIlVl 




54157 


.0 


147 


,2670 ± 


0. 


1061 


TXT 1 ' 1 VI AT /TC/^FJT 

IN IbuKAb/lSCjRl 


25 


54160 


.0 


147 


.2800 


+ 





.0755 


TNTFGRAT /TSGRT 25 

1.1 1 L L j V I IV I V L j/ 1i.JV.I1V1 




54163 


.0 


147 


,3790 ± 


0. 


0746 


TXT 1 ' 1 VI AT /11VD1 

IN IbuKAb/lSuKl 


25 


54166 


o 


147 


.4210 


+ 


o 


0999 


TNTFGRAT /TSGRT 25 

1.1 1 L I j V 1 1 V I V 1 it 1 v. I V.I 




54169 


.0 


147 


,4370 ± 


0. 


0758 


TXT 1 ' 1 .' VI AT / 1 1 • /" ' T~l I 

IN IbuKAb/lSCjRl 


25 


54175 


o 


147 


.5010 


+ 





.0742 


TNTFGRAT /TSGRT 25 

111 L L j V I IV I V 1 -/ IklVIlVl 




54178 


.0 


147 


,5084 ± 


0. 


0693 


TXT 1 ' 1 VI AT /I tv D I 

IN IbuKAb/lbuKl 


25 


54187 


o 


147 


.6360 


+ 


0. 


.0680 


TNTFGRAT /TSGRT 25 

1.1 1 I I j VI 1 V I V Lj/ IklVIlVl 




54193 


.0 


147 


,7020 ± 


0. 


1296 


TXT 1 ' 1 VI AT /TC/^FJT 

IN IbuKAb/lSCjRl 


25 


54196 


.0 


147 


.7520 


+ 





.0915 


TNTFGRAT /TSGRT 25 

jLL 1 L I j V 1 1 V I V 1 -/ IklVIlVl 




54205 


.0 


147 


,6560 ± 


0. 


5624 


TXT 1 ' 1 VI AT /11VD1 

IN IbuKAb/lbuKl 


25 


54208 


o 


148 


.0240 


+ 


0. 


1515 


TNTFGRAT /TSGRT 25 

L.L 1 I IjVllVy VLj/ IklVIlVl 




54211 


.0 


147 


,8410 ± 


0. 


1414 


TXT 1 ' 1 .' VI AT / 1 1 • /" ' F) I 

IN IbuKAb/lSCjRl 


25 


54214 


.0 


147 


.9820 


+ 





.0789 


TNTFGRAT /TSGRT 25 

1.1 V L L j V I IV I V 1 -/ IklVIlVl 




54334 


.0 


149 


1250 ± 


0. 


1883 


TXT 1 ' 1 VI AT /I tv D I 

IN IbuKAb/lbuKl 


25 


54337 


o 


149 


.2490 


+ 


0. 


1609 


TNTFGRAT /TSGRT 25 

1.1 1 I I j V 1 1 V I V Lj/ IklVIlVl 




54340 


.0 


149 


,2190 ± 


0. 


1271 


TXT 1 ' 1 VI AT / 1 1 • /" ' ri I 

IN IbuKAb/lSCjRl 


25 


54346 


.0 


149 


.2859 


+ 





1262 


TNTFGRAT /TSGRT 25 

111 L 1 j V 1 1 V I V 1 -/ IklVIlVl 




54355 


.0 


149 


,4520 ± 


0. 


1408 


TXT 1 ' 1 VI AT /TC/^FJT 

IN IbuKAb/lSuKl 


25 


54358 


o 


149 


.2110 


+ 





1932 


TNTFGRAT /TSGRT 25 

1.1 1 I IjVllVy VLj/ IklVIlVl 




54361 


.0 


149 


,7855 ± 


0. 


4506 


TXT 1 ' 1 VI AT /TC/^FJT 

IN IbuKAb/lSCjRl 


25 


54364 


.0 


149 


.5520 


+ 


o 


.4772 


TNTFGRAT /TSGRT 25 

XL 1 L L j V I IV I V 1 -/ IklVIlVl 




54367 


.0 


149 


,7460 ± 


0. 


2236 


TXT 1 ' 1 VI AT /TC/^TUT 

IN IbuKAb/lbuKl 


25 


54370 


o 


149 


.3710 


+ 


0. 


1449 


TNTFGRAT /TSGRT 25 

ll'l I I j V 1 1 V I V Lj/ IklVIlVl 




54379 


.0 


149 


,6473 ± 


0. 


1996 


UN 1 J_iVJlVrVJ_,/ 1 k I v J Iv 1 


25 


54388 


.0 


149 


.7940 


+ 





1169 


TNTFGRAT /TSGRT 25 

ll'l L I j V 1 1 V I V I _// IklVIlVl 




54391 


.0 


149 


,6390 ± 


0. 


1620 


INTEGRAL/ISGRI 


25 


54541 


o 


150 


5900 


+ 


o 


.3070 


TNTFGRAT /TSGRT 25 

ll'l I IjVllVy VLj/ IklVIlVl 




54704 


2 


152 


,7987 ± 


0. 


0017 


Fermi/GBM 25 




54719 


5 






+ 


o 


.0027 


Fermi/GRM 25 

X ^.11111/ VJ1J1V1 




54728 


1 


l 


Ulj^ ± 


0. 


0037 


Fermi/GBM 25 




54736 


1 


153 


0986 


+ 





.0024 


Fermi/GBM 25 




54744 


.1 


153 


1724 ± 


0. 


0018 


Fermi/GBM 25 




54752 


1 


153 


.2497 


+ 





0031 


Fermi/GBM 25 




54769 


.0 


153 


,3928 ± 


0. 


0030 


Fermi/GBM 25 




54775 


,9 


153 


.4629 


+ 





0016 


Fermi/GBM 25 




54784 


.0 


153 


,5438 ± 


0. 


0012 


Fermi/GBM 25 




54791 


,9 


153 


.6298 


+ 





.0012 


Fermi/GBM 25 




54800 


.4 


153 


,7102 ± 


0. 


0024 


Fermi/GBM 25 




54808 


1 


153 


.7824 


± 


0. 


.0025 


Fermi/GBM 25 




54815 


.7 


153 


,8466 ± 


0. 


0021 


Fermi/GBM 25 




54824 


.5 


153 


.9263 


+ 





0010 


Fermi/GBM 25 




54831 


.9 


154 


,0113 ± 


0. 


0007 


Fermi/GBM 25 




54840 


,0 


154 


1071 


+ 





.0007 


Fermi/GBM 25 




54847 


.9 


154 


1913 ± 


0. 


0011 


Fermi/GBM 25 




54856 


.2 


154 


.2690 


+ 





.0022 


Fermi/GBM 25 




54880 


.0 


154 


,4594 ± 


0. 


0030 


Fermi/GBM 25 




54887 


,9 


154 


.5198 


± 


0. 


.0023 


Fermi/GBM 25 




54904 


.7 


154. 


.6587 ± 


0. 


0032 


Fermi/GBM 25 




54912 


.2 


154 


.7181 


+ 





.0020 


Fermi/GBM 25 




54919 


.8 


154 


,7824 ± 


0. 


0020 


Fermi/GBM 25 




54927 


,9 


154 


.8590 


+ 





.0024 


Fermi/GBM 25 




54936 


.3 


154 


,9303 ± 


0. 


0016 


Fermi/GBM 25 




54944 


.0 


154 


9999 


+ 





0016 


Fermi/GBM 25 




54952 


.1 


155 


,0761 ± 


0. 


0021 


Fermi/GBM 25 
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Table B.l. Pulse period measurements of GX 1+4 (continued). 



MJD 




1 spiny'* J 




Instrument 


MJD 








Instrument 




54959 


5 


155 1347 + 


0.0023 


Fermi/GRM 2i 

A 1111/ V J 1 J • V 1 


54968 


o 


1 55 2070 + 


0.0019 


Fermi/GRM 

A ^ 1 1111/ VJ1J1V1 


25 


54975 


9 


1 55 2761 + 


0.0021 


Fermi/GRM 25 

A Li 1 1 11/ VJ1J1V1 


54984 


o 


1 1446 + 


0.0014 


Fermi/GRM 

A L 1 1 111/ v 1 1 J 1 » 1 


25 


54992 


o 


155.4194 ± 


0.0010 


Fermi/GRM 25 

X Islllll/ V 1 1 J 1 V 1 


55000 


o 


1 4976 + 


0.0010 


Fermi/GRM 

A L L 1 1 11/ VJ1J1V1 


25 


55007 


9 


1 55 5781 + 

1 _ ' *J . —/ 1 O 1 _!_ 


0.0010 


Fermi/GRM 25 

X 1^11111/ V J J J 1 V 1 


55016 


1 


1 55 6549 + 


0.0019 


Fermi/GRM 

A w 1 1111/ VJ1J1V1 


25 


55023 


9 


155.7235 ± 


0.0018 


Fermi/GRM 25 


55031 


9 


155.7919 ± 


0.0023 


Fermi/GRM 


25 


55040 





155.8628 ± 


0.0026 


Fermi/GRM 25 


55047 


9 


155.9249 ± 


0.0016 


Fermi/GRM 


25 


55056 





156.0042 ± 


0.0009 


Fermi/GRM 25 


55063 


8 


156.0889 ± 


0.0007 


Fermi/GRM 


25 


55072 


3 


156.1731 ± 


0.0015 


Fermi/GRM 25 


55079 


8 


156.2443 ± 


0.0016 


Fermi/GRM 


25 


55087 


9 


156.3116 ± 


0.0021 


Fermi/GRM 25 


55120 





156.5714 ± 


0.0011 


Fermi/GRM 


25 


55128 


1 


156.6515 ± 


0.0007 


Fermi/GRM 25 


55135 


9 


156.7381 + 


0.0008 


Fermi/GRM 


25 


55143 


9 


156.8283 ± 


0.0007 


Fermi/GRM 25 


55151 


9 


156.9130 ± 


0.0009 


Fermi/GRM 


25 


55159 


9 


157.0009 ± 


0.0006 


Fermi/GRM 25 


55168 





157.0931 ± 


0.0005 


Fermi/GRM 


25 


55176 





157.1903 ± 


0.0005 


Fermi/GRM 25 


55184 


2 


157.2875 + 


0.0006 


Fermi/GRM 


25 


55191 


9 


157.3714 ± 


0.0009 


Fermi/GRM 25 


55200 





157.4465 ± 


0.0014 


Fermi/GRM 


25 


55207 


9 


157.5144 + 


0.0022 


Fermi/GRM 25 


55224 





157.6355 + 


0.0023 


Fermi/GRM 


25 


55232 





157.7167 ± 


0.0006 


Fermi/GRM 25 


55237 


5 


157.7838 ± 


0.0019 


Fermi/GRM 


25 
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